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ABSTRACT 


A numerical study has been made to develop a technique for 
studying displacement body and curvature effects on incompressible and 
compressible laminar bo^lndary layers. 

This technique has been applied to study such effects on incompres- 
sible flow around cylinders at moderate to low Reynolds numbers and for 
compression ramps at hypersonic Mach numbers by employing a finite-dif- 
ference method to obtain nmerical solutions. 

The results indicate the technique can be applied successfully in 
both regimes and does predict the correct trend in regions of large 
curvature and displacement body effects. 

This study concludes that curvat-ure corrections should only be 
attempted in cases where all displacement effects can be fully 


accounted for. 
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I. FORMULATION OF THE PROBLEM 


1. 1 Introduction 

Since the 1950' s, much attention has been given to higher-order 
boundary- layer theory due to the academic challenge it presents and 
also in hope that this inclusion of higher -ordered terms would enable 
one to extend the boundary-layer concept to flow problems which 
normally would be considered to exceed the limits of classical or 
first-order theory. 

Even though higher -ordered terms were beginning to be calculated 
in the 1950' s, it was only recently, in the 1960's, that all second- 
order contributions were correctly identified. This was accomplished 
by Van Dyke (1962a, 1962b) for incompressible theory and independently 
by Van Dyke (1962c), Maslen (1965)> and Lenard (1962) for compressible 
theory, all of whom employed the method of matched asymptotic expan- 
sions which results in two separate sets of governing equations - one 
for the classical or first -order contribution and another for the 
second -order contributions. 

Effects to be studied here which were found to be of higher order 
in an asymptotic sense can generally be classified as arising either 
due to the surface curvature of the body under consideration or due 
to its boundary layer interacting with the external flow field. Of ' 
the two second -order effects which will appear in this study, perhaps 
the easiest to understand and usually the most difficult to calculate 


1 
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is what is termed the displacement-body or displacement- speed effect 
which arises due to the boundary layer interacting with the external 
mainstream. This effect appears as a result of the viscosity influence 
near a solid surface which tends to slow down the fluid in the region 
near the wall and force it outward so that the effective body presented 
to the mainstream is the original body plus the displacement thickness 
of the boundary layer. Since the boundary- layer thickness is not 
known beforehand, the difficulty in computing this effect is immediately 
recognized. The other higher-order effect considered here is that of . 
surface curvat\ire and presents no unusual difficulties. 

In addition to clearly indicating what the second -order contribu- 
tions are, another advantage of filtering out what are normally higher- 
order effects as done by Van Dyke, Maslen, and Lenard is that all higher- 
order equations are linear and therefore can be divided into several 
additive effects each of which has a simpler physical meaning. However, 
the purpose of this study is to examine problems where there are large 
displacement and curvature effects. Thus, the distinction between a 
first-order and a normally second-order effect may be small and as such 
some advantage may be gained if we do not try to separate them as 
normally done but rather seek only to filter out those strictly third - 
order and higher contributions which appear in the analysis. The final 
equations will therefore contain both the first-order and second -order 
contributions given by Van Dyke, Maslen, and Lenard plus some addi- 


tional third -order terras. 
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Most studies of higher-order effects have been made adopting the 
equations which have the higher-order contributions filtered out. In 
reviewing the literature, we examine subsonic and supersonic external 
flows separately since there appear basic differences in the nature 
of their governing equations which affect the ease with which displace- 
ment effects can be computed. Attention here is paid to cases where 
displacement effects have been computed since once these are known 
the effects of surface curvature can be easily handled. For subsonic 
theory, the effect of displacement speed cannot usually be calculated 
directly owing to the fact that the external governing equations are 
elliptic thereby making its effect dependent upon the entire course of 
the boundary layer. Thus, there have been very few cases where dis- 
placement effects have been computed in subsonic theory. Kuo (1953) 
calculated the displacement effect for flow over a finite flat plate. 
Few shapes other than the finite flat plate have been investigated, 
of these Van Dyke (1964) treated the parabolic cylinder where he was 
able to calculate the displacement effect only near the leading edge. 
Devan (1964) treated the Rankine heilf-body numerically. 

In supersonic or hypersonic flow, the inviscid or outer flow 
equations are hyperbolic thereby making the displacement effect local 
in influence and therefore simplier in theory to calculate. Maslen 
(1952) treated the semi-infinite plate with its mixed subsonic and 
supersonic external flow regimes employing linearized subsonic and 
supersonic inviscid theory for the flow due to displacement effects. 
Most investigations of compressible second-order effects have been 
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applied to blunt bodies in hypersonic flow. Individual second- order 
effects at a stagnation point have been calculated by Lenard (1962)^ 
Davis and Fliigge-Lotz (1964a), and Fannelop and Fliigge-Lotz (1965), 
whereas extended solutions around blunt bodies were first reported by 
Davis and Fliigge-Lotz (1964b) for a sphere, paraboloid, and hyperboloid, 
and by Fannelop and Fliigge-Lotz (1964) for a circular cylinder. 

Of particular interest in this study are bodies with sharp leading 
edges which are placed in a hypersonic mainstream. For problems of 
this nature, there arise very large displacement effects - particularly 
near the leading edge. Hayes and Probstein (l959) reviewed the early 
literature relating to this subject which, due to the strong viscous 
interaction present, is usually treated separately from supersonic 
flows. There are many approximate solutions for flat plates in hyper- 
sonic flow appearing in the literature. Dewey (1965) adopted the local 
similarity concept for his analysis. Later, Chan (1966) employed the 
momentum-integral method in his study. More recently, Sullivan (1969) 
employed Lees (1956) cold-wall similarity analysis to extract an 
approximate solution. Only one exact solution^ has been reported for 
a flat plate in hypersonic flow and that was by Fliigge-Lotz and 
Blottner (I962). Only in one case, has this theory been extended to 
cover bodies with surface curvature and this was reported by Stollery 
(1970) who did not include surface curvature effects in his approxi- 
mate analysis. Thus to the author's knowledge, no studies have been 

^Exact in the sense that no approximations were made to the 
governing bound ary- layer equations. 
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made in this area which include both displacement and curvature 
effects. 

The first problem to be investigated in this study is that of 
flow around a circular cylinder at moderate to low Reynolds numbers. 

This flow will be investigated employing the composite second-order 
incompressible boundary-layer equations and it is felt that displace- 
ment and curvature effects shoiald be sizable for the Reynolds number 
cases studied. The second flow problem to be studied is that of flow 
up a compression ramp in hypersonic flow. It has been noted by Van Dyke 
( 1969 ) that even for high Reynolds numbers, the effects of displacement 
and surface curvature become significant for large Mach number flows. 

For this study, the compressible second-order boundary-layer equations 
will be employed. 

From this study, large displacement and curvature effects were 
found to exist. However, when both displacement and curvature effects 
were accounted for, the differences observed in the wall properties 
from that predicted by classical boundary- layer theory were only 
noticeable in regions where the displacement effect could no longer 
be classified as a higher-order effect and this occurred near the 
leading edge of the compression ramp in hypersonic flow. 

1 . 2 Formulation of the Governing Equations 

For flow of a viscous fluid past a stationary body, Prandtl (1904) 
observed that as the Reynolds number becomes large, a thin layer 
develops' near the body surface where the effects of viscosity and 
heat transfer are significant and outside of which they are negligible. 
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> 

Prand'tl formalized this concept from an order of magnitude analysis 

of the full Navier- stokes equations and thereby produced the now 

i 

famous boundary-layer equations. Following this lead, Van Dyke (1962a) 
used the method of matched asymptotic expansions to obtain a formal 
statement of the boundary-layer concept valid to second order in the 
inverse square root of Reynolds number. Analyses of the individual 
first and second- order effects for several incompressible flows were 
made by Van Dyke (l962a) and (1962b). Siniiliar analyses for compres- 
sible flow were made independently by Maslen (1965)^ Lenard (1962), 
and Van Dyke (1962c). Applications of these are discussed in 
Van Dyke (1969)* Here we choose not to filter out all higher-ordered 
terms or to separate the second-order contribution from the first 
order Iraii.dtl boundary-layer equations, iiie purpose of this approach 
is to determine a composite set of equations which formally will 
contain both the .Vraiidt.l boundary-layer equations and the second-order 
contributions given by Maslen (19(>3^ » Denard (1962), and Van Dyke 
(1962c 'I . 

The analysis'*' begins by defining the perturbation parameter 

e - 1/ V Kc , (1.2-la) 

and thereafter stretching the normal coordinate and velocities 
according to the relation 

'''The following analysis was first presented by Davis, Werle, and 
Wornom (1970) and is presented here as a review and in order to show 
m.ore of the details. 
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n = eN, 


(1.2-lb) 


and 

V = ev, (l.2-lc) 

with all other variable ass\:iined to be order one. 

These are now introduced into the Wavier -Stokes equations given 
by Van Dyke (1962c) with all terms of 0(l) and 0(e) kept so that 
the governing equations valid to second order (i.e., 0(e)) become: 


CONTINUITY EQUATION 


(r-^pu)g + (hrJpv)j^ = 0 (exact) 


( 1 . 2 - 2 ) 


where d = 0 for plane flow and j = 1 for axisymmetric flow; 


S-MOMENTUM EQUATION 


u 


Uo + v(hu). 


6s hr'^ oW 


(l.2-3a) 


where the shear stress is given by 


T = |j.(ujj - h]\]u/h) 


(l.2-3b) 


N-MOMENTUM EQUATION 


^N 2 
p = ^ pu 

N h 


(1.2-4) 
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ENERGY EQUATION 


p(iiHg + hvHjj) 


1 8 
8n 


hr^ ( q + UT ) 


where 


q = tL ^ j (heat transfer) 
a dN 


H = T + i , (total enthalpy) 


and 


P*C* 

a - . ( Prandtl number ) 

k* 


In the above equations the coordinate scale factors r and 
defined as 


r 


+ eN cos 0 


wall ’ 


and 


h = 1 + gkN . 


The viscosity law employed here will be either 
a. Sutherland's law 


d 


= 1 + C 

T + C 


where 


C* 

(7 - 1) M^T* 
' 00 00 


(l.2-5a) 

(l.2-5b) 

(1.2-5C) 

(l.2-5d) 
h are 

(l. 2-6a) 
(l.2-6b) 

(l.2-7a) 


C 




(l.2-7b) 
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and C* = 198. 6°R for air, or the 
h. Linear viscosity law 

JL = c(s)-^-, (l.2;-8a) 

M-OO *^00 


where |i^ is calculated using the Sutherland viscosity law and C(s) 
is determined by requiring the linear law to match Sutherland's law 
at the wall, that is. 


C(s) 


1 + ^ 
T^/T^ + C*/T* V ‘ 


(l.2-8b) 


Only cases of constant wall ten^erature will be studied here so that 
here 

C(s) = Coo = constant. (l.2-8c) 


The boundary conditions represent the no- slip condition on 
velocity 

u(s,o) =0, (1.2. -9a) 

the injection velocity at the boiindary 

v(s,o) = v^(s) , (l.2-9b) 

and the total enthalpy at the wall 


H(s,o) = H^(s) , 


(l.2-9c) 
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or the adaibatic wall condition 

/IS'l =0. (l.2-9d) 

V an/wall 

The outer boundary conditions are more subtle than the wall conditions 
because the edge of the boundary layer is not well defined. For this 
reason we must look for matching conditions to mate the viscous flow 
to the outer inviscid flow state. 

In the present investigation, only cases where the outer flow 
was irrotational were studied. Thus the first matching condition will 
require that as N goes to infinity (i. e. , the region where the 
boundary layer merges with the inviscid flow), the boundary-layer 
vorticity should asymptotically approach the inviscid value. 
Mathematically stated 

CO ~ 0 , as N 00 , (l. 2-10 ) 

where to second order the vorticity in the boundary layer is given by 


CD = i h ^(hu) 


N ’ 


( 1 . 2 - 11 ) 


"^The fact that the outer boundary condition must be applied at 
infinity can be seen by letting n^ be a point in the region where 
the boundary layer merges with the inviscid flow. Then the matching 
conditions are applied as 


N Nq = no/e . (l.2-10a) 

Since boundary-layer theory becomes exact in the limit as e ^ 0, 
then N ^ 00 at the outer reaches of the viscous region. 
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From equations (1.2-10) and (1.2-11), we find that 

(hu)|g 0 as N 00 , (1.2-12) 

which results in 

hu ~ f(s) as N ^ oo , (1.2-13) 

To determine the function f(s), we apply the matching condition that 
the viscous u-velocity profile as N becomes large matches the 
inviscid profile near the body surface. This is the same matching 
condition as used by Murphy (1953)^ Van Dyke (1969), and Davis, 
Whitehead, and Wornom (1970). This implies that 

f(s)/h ~ U(s, n) for small n, (1.2-l^f) 

where U(s,n) is the inviscid U-component of velocity. Noting that 
h = 1 + ten , we cau expand both sides of equation (l.2-l4) in a 
Maclaurin series for small n to obtain 

f(s)(l - Kn + ...) ~U(s,o) n + ... . (1.2-15) 

y onjs, o 

From the condition that the outer inviscid flow be irrotational comes 
the relation that 

S(hU ) BV 
3n 5s 


(1.2-16) 
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Now if we solve equation (l.2-l6) for ^ and substitute it into 
equation (1.2-15), we find that 


SV\, 


f (s) (l - Kcn + . . . ) ~ U(s,o) (l - Kn) + + . . . . for small n . 

(1.2-17) 

In the outer flow the V-component of velocity can be expanded as 

+ 




V = + 6V2 + . . . 


(1.2-18) 


Thus equation (1.2-17) can be expressed as 


f(s)(l 


Kn + ...) ~ U(s,o)(l - Kn) + 


/v» <v> 



(1.2-19) 


Noting that U(s,o) = Ug(s), ( 8 Vi/ 8 s)g^o = 0 for any body, and 
n = 0(e) for proper matching, we find from equation (1.2-19) that 
to second order 


f(s) = Ue(s) . (1.2-20) 

Here Ug is the inviscid s-urface speed on the body under considera- 
tion and must be valid to second order, that is, it must contain the 
displacement- speed effects delineated by Van Dyke (1962c). Therefore, 
the resulting matching condition on the u-con5)onent of velocity 
becomes 

"'"See Van Dyke (1962c). 
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u ~ Um = Ug/h as N -> “ • (1.2-21) 

The outer boundary condition on enthalpy can be found as follows. 

As N ->oo in equations (1.2-3) and (l.2-5)> it will be required that 
the viscous terms vanish (i.e., the governing equations yield the 
inviscid flow equations ) . This gives the condition that 

p(Ujj^Hg + hvHjj) ~ 0 as N -> oo . (1.2-22) 

Letting H be denoted as as N -> «> , it is seen that 

Hjn = constant (1.2-23) 

is a solution to equation (1.2-22). The obvious choice for the con- 
stant is 

= Hg, (1.2-24) 

where Hg is the inviscid total enthalpy on the body surface here 
taken to be constant. 

Therefore, the matching condition on total enthalpy becomes 

H ~ = He as N oo . 


(1.2-25) 
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In addition to the boundary and matching conditions given by equa- 
tions (1.2-9), (1.2-'21), and (1.2-25), we must also identify how the 
following quantities match with the inviscid flow. As N “ the 
viscous terms must approach the following inviscid matching profiles, 
denoted by the m subscript, 

q ~ Pm> (l.2-26a) 

T~Tjn, (l. 2 - 26 b) 

P ~ Pjn> (1.2-26 c) 

and 

p~p^. (l.2-26 d) 

Letting N become large in equation (l. 2-5c), we find after substitu- 
tion from equations (1.2-21) and (1.2-25) that 

T„~He-u^/2. (1-2-27) 


Similarly, equations (1.2-7) and (1.2-8) give that 




1 + C 
Tm + C' 


T, 


3/2 


m 


(l.2-28a) 


or 


^‘•m _ fi ■^ra 

— ^OO • 


(l.2-28b) 
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Letting N approach infinity in the s-raomentum equation and noting 
that the viscous terms will he required to vanish as N <» , we obtain 


Pra^m-§^ ^ ^ asN^oo. 


(1.2-29) 


If we now substitute 



(l.2-29)j we find that 


from equation (1.2-27) into equation 


^ = a 

as ^ as 


(1.2-30) 


Making the substitution for in equation (1.2-30) from the equation 

of state, we find upon integration that 


Pm(s^N) 


Ti„(s,N) 


- Z ,. 

7-1 

W(N) 


as N -> 00 . 


(1.2-31) 


To determine the function W(W), we apply the matching condition that 
the viscous press\ure distribution for large N should match the 
inviscid pressure for small n, that is, 


Pjjj(s,N) ~ P(s,n) for small n, (1.2-32) 

where P(s,n) is the external inviscid pressure. 

For small n we can express P(s,n) in a Maclaurin series about 


n = 0. Thus 



l6 


I^(s,n) 


P(s,o) + 


3p 


n + 


s, o 


From Van Dyke (1962c), the n-momentum equation evaluated at 
becomes 


Sp'- 


= -Pe 


U. 


s, o 


6v 


s,o 





where Ug = U(s,o), Vg = V(s,o), and Pg = p(s,o) . This can 
be substituted into equation (1.2-33) to obtain 


P(s,n) = Pg - pg 



+ V, 


3n 


- KUg 


n + 


s, o 


where 


Pe(s) = P(s,o) . 

To second order within the boundary -layer region, equation ( 
may be written as (using the equation of state) 


P(s,n) 



Noting that 


Tm = (Tg + Ug/2) - u^/2 , 


(1.2-33) 
n - 0 

(1.2-3^) 

now be 

(l.2-35a) 

(l.2-35b) 

.2-35) 

(1.2-36) 


(1.2-37) 
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equation (1.2-31) can be written as 


= W(N) 




h2 



which when expanded for small n gives 


(1.2-38) 


p„ = W(N)T, 


7 

7-1 


1 

7 - 1 Tp 


Kn 



(1.2-39) 


If we now apply the matching condition that pjn(s,N) ~ P(s, n) for 
small n, we find that equation ( 1 . 2 - 36 ) and (1.2-39) give that 


7 

W(N) = W = Pg/T^"^ = constant (l.2-U0) 


for irrotational flow. Thus, equation (1.2-31) be written as 


P^(s,N) 



(1.2-41) 


where Pg and Tg are the inviscid pressure and temperature on the 
body under consideration. Likewise, from the equation of state, 
p^(s,N) can be written as 

pJs,W) = , (1.2-42) 

where Pg is the inviscid density on the body under consideration. 
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It was noted previously that the viscous terms would be required 
to vanish from the energy and s-momentum equations as W oo . As 
the equations now stand this is not true. To see this, let N 
approach infinity in the s-momentum equation, where the right-hard 
side (RHS) becomes 


RHS ~ 


1 _ 8 _ 
hrj S'* 



r\ 


h h / 


as N oo ^ 


(1.2-43) 


which when expanded gives 


RHS 


hrJ 






(1.2-44) 


Noting that h = 1 + €kN, r = r^ + eN cos ©wall’ '^N “ ^N = 

we see that although the viscous terms do not vanish, they can be 
seen from equation (1.2-44) to go to third order. Since the governing 
equations are valid formally only to second order, these third- order 
terms which are bothersome to the numerical analysis can be corrected 
for by subtracting the limiting form of these same terms from the RHS. 
Applying this correction the s-momentum equation becomes 


P I^U Ug + V 


(hu) J 


8p _ 1 8 

8s hrj 8N 


2 i / 
h r'Jq{Uj^ 


% u /i ^ 

h I ~ u njj ’ 


(1-2-45) 
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Similarly, letting N approach infinity in the energy equation we 
find that the RHS becomes 


RHS 


rJ 


Sn 


ih 


hr ^ (1 - 2o) 

_ a 1 


as W ->oo, (l.2-46a) 


which after expansion becomes 


RHS ~ — 





.2 

+ r3 — 


dp 


dT 


u 

m 


m 


2rJ 


p V? 

m m 


I]- 

(1.2-U6b) 


Again, the RHS does not vanish as required but it can be seen from 
equation (l.2-46b) to be third order and as such will be corrected 
for the same manner as the s-moraentum equation. Applying this correc- 
tion, the energy equation becomes 


p(uHg + hvHjj) 


= i A [ ~hr’^p ,! ^dH _( 


•(1 - a) u ^ 
pJ dN ^ a l^dN dN 


cr u 


2 % 
h 


(1 - 2o)^^ u^l 

p h "y 


(1.2-47) 


One final modification which leads to a simplier computational 
scheme involves writing the pressure gradients in equations (1.2-4) 
and (1.2-45) in terms of the velocity profile by first integrating 
equation (1.2-4) from any point Nq to an arbitrary point W to give 
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p(s,N) 


p(s,Wo) 



dW . 


Realizing that 


p™(s,N) = p (s,No) + 


N 


hi 




.2 


N, 


~ dW , 


o 


equation (1.2-48) can be written as 


(1.2-48) 


(1.2-49) 


p(s,W) 


Pra(s>N) 


+ 


p(s,Nq) 


Pm(s^No) 


W, 


N 


° Mp 
h I m m 



dN . 


(1.2-50) 

n_ 

In the limit as e -> 0 for n^ fixed, Nq = -x -> oo and 
p(s,Nq) -^Pj^(s,Nq) , so that equation (1.2-50) becomes 


p(s,W) 


Pm(s>N) 



i^/p u2 
h m m 



(1.2-51) 


where p^(s,W) is given by equation (l.2-4l). 

1. 5 Summary of Governing Equations in Boundary- Layer Variables 
The compressible boundary-layer equations valid formally to 
second order are 
the CONTINUITY EQUATION 



(1.3-1) 
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the S -MOMENTUM EQUATION 


Su. 


3[u Us + v(hu)]^] - PjnU^ 


1 a 


j 


h^rJJuN - u (l - 




N 


( P u 

h I m m 


pu" 


2^ra 


dN 


(1.3-2) 


and the ENERGY EQUATION 


p(\iHs + hvHN) = 


^ A 

r I SNi 


hr Jp r^H / 

a ^Sn 


N 

- o)u — - a u 

Sn 


-(1 - 2 a ) 


tiiiiusl 


M. h 


■m; 

J 


2 

h 


(1.5-5) 


The pressure in the boundary layer is given by 


p(s,N) = p^(s,N) - pu^jdN (l.5-i+) 


h 


EQUATION OF STATE 


p(s,N) = 2 LjlJ. p(h - u^/2) 
7 


(1.3-5) 


VISCOSITY LAWS 

Two viscosity laws will be employed. They are: 

a. The Sutherland viscosity law 

,, _ 1 + C .nV2 


T + C 


(l.3-6a) 
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where 


c* 

c 1 — ^ 

(7 - 1 ) ' 

(l.5-6b) 

and 


b. The linear viscosity law 


U = 

(l. 3 - 7 a) 

Too 


where 


^ _ 1 + c */ t * . 

“ T^/T^ + C*/T^ V T^ 

(l- 3 - 7 b) 

In equations (l. 3 - 6 b) and (l. 5 - 7 b), C* is taken 

to be 198. 6°R 

for air. 


BOUNDARY CONDITIONS 


The boundary conditions are the no- slip condition 


u(s, o' = 0 , 

(1. 5-^a' 

the surface injection condition 


V(s, 0 )= V^(s) , 

(l.3-8b) 

a specified wall temperature 


H(s,o) = H^^(s) , 

(1.3- 8c) 

or the adaibatic wall condition 


( 3 h/ 8 n)^^ = 0. 

(l.3-8d) 


The matching profiles which constitute the outer boundary conditions 
are 

Uj^ = Ug/h as N -> 00 ^ 


u(s,W) 




(l. 3 - 8 e) 
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H(s,W) ~ Hjn = Hg as N , (l.3-8f) 

_2L_ 

‘V — 1 

p(s,N) ~ pjn = Pe(T^/Te) as N ~ , (l.3-8g) 

1 

p(s,N) ~ Pm = Pe('^m/'I’e)^ ^ as W oo , (X. 3-8h) 


and 

T(s,N) ~ = Hg - Uj^/2 as N ^ . (l.3-8i) 

In the above equations Hjjj, Pm> Pm^ inviscid 

profiles that the viscous profiles must match as the boundary layer 
merges vdth the inviscid flow region. 

1. 4 The Governing Equations in Similarity Variables 

It will prove convenient later if the governing equations are 
put into similarity form so that when a self-similar solution exists 
the equations will reduce to the proper form. The similarity form 
is also better suited for niimerical solution since the boundary layer 
does not grow as fast in these coordinates. 

Following the example of Hayes and Probstein (1959) ^ Blottner 
(1970), and Davis, Whitehead, and Wornom (197O), we choose the trans- 
formed independent coordinates 5 and t) to be 




2h 


(1.4-1) 
(1.4-2) 

Equations (l. 4-1-2) are generalized forms of the Gortler and Howarth- 
Dorodnitsyn transformations. Note that for the incompressible case 
p = Pg = |ig = 1 and equations (l. 4-1-2) reduce to the expressions 
employed by Davis, Whitehead, and Wornom (1970) and for classical or 
first-order theory h = 1 and r = Tq so that equations (1.4-1) and 
(1.4-2) reduce to those coordinates originally due to Levy (1954) and 
Lees ( 1956 ) and adopted by Blottner (1970), and others. 

The transformed tangential component of velocity is defined as 

^ = u/um ^ (1-4-5) 

and the transformed total enthalpy is defined by 

G = H/Hm. (1.4-4) 

Applying the above transformations to the continuity equation we 
obtain upon integration that 

hr^pv = Bs /gT J F dqj - |2 \| 2 ^ F + Pw^q^v^ . 



(1.4-5) 
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The transformed v- component of velocity will be chosen to be 


V = - \f2i A j F dT) j + p^r8v^\/2t/(d|/ds) . 

(l. 4-6a) 

Thus, in terms of the similarity V-component of velocity the physical 
v-component becomes 


(l.4-6b) 


= -K- ( — — ^ ^ \/ 2| F ) . 

pr^h V ds y 2? ds 


Differentiating equation (l.4-6a) the transformed continuity equation 
can be written as 


A+2|~+F=0 
dr] d| 


(l. h-Gc ) 


In the new coordinate system the convective operator becomes 




dN, 


\ _ PUm d| 

25F A + V A 

/ 25 ds 

d5 dq 


(1.4-7) 


Applying the above transformations to equations (l. 3-2-3) the trans- 
formed second -order compressible boundary-layer equations are obtained 
as 

CONTINUITY EQUATION 


A + 2| — + F= 0 
dq d^ 




(1.4-8) 
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S -MOMENTUM EQUATION 




St] ] 


6t] 




F - 


h \ 


P^m 

Pm^ /-I 


h" ^ - — ^ + 2|F — + 

\ P / 


.V|.f2p^ .2t|l,/p, 


/ 


(l.i^-9) 


ENERGY EQUATION 


A 

6t) 


1 He 


2r 

(a 


1)F M + (1 _ 2a) ^ 

6t) h V p^Z 


I - V Ij . 


(l.i^-lO) 


The pressure at any point in the boundary layer is given by 


p(s,N) = p„ + u^ I 


(1.4-11) 


where u^^^ and p^ are given by equations (l. 5-8e) and (l. 5-8g). 
EQUATION OF STATE 


p = 21^ p(HgG - u2/2) . 


(1.4-12) 


Quantities not defined in equations (1.4-8-12) are defined as 


R = r/ro , 


(1.4-13) 


I = (p^y (PgPe) » 


(1.4-14) 
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and 


21 aufli _ 2g dUg 2g 5h 

Ujn Ug d| h dl, 



pr 




Pm / 


■ pil :> 


= (Prat^m)/(Pe^^e) ' 


Equations (l. 3-6-8) supply the remaining relations needed 
a solution. 

BOUNDARY CONDITIONS 

The boundary conditions are the no slip condition 

F(5,o) = 0, 


a specified wall temperatirre 


G(l,o) = 


or zero heat transfer 




the surface blowing velocity 


(iA-13) 

(1.4-16) 

(1.4-17) 
for computing 

(l.4-l8a) 

(l.4-l8b) 

(l.4-l8c) 


VU,o) = Vji) , 


(l.4-i8d) 
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the outer matching condition on velocity 

F~1 as N->oo, (l.4-l8e) 

and the outer matching condition on total enthalpy 

G ~ 1 as N ^ 00 . (l.4-l8f) 

!• 5 Bound ary- Layer Properties 

The dimensionless coefficient of skin friction is defined as 


''^wall / Sf\ ^ ^ ^ 

= T = ^ ^ (1.5-la) 

|/2I V^Vwall 


where 


r wall 


wall 


(1.5-lb) 


This leads to the usual definition of skin friction 


Cf = Cf/e = 

)/ 2? \^’^/wall 


(1.5-2) 


The dimensionless 
Stanton number defined 


coefficient of heat transfer is given by the 
as 


q'^wall _ eHeP^^rpJ'Ue / (1.5.3') 

P»U*(H* - H*) a]jW {To - Tw)\^^4all' 
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where 


q* 

^ wall 


^ ^>^*/wall ’ 


(1.5-4) 


We can remove the Reynolds niomber dependence from the Stanton niomher 
by defining 


St^ = Stoo/e 


HePw*"w^o^Ue / ba\ 

^ (Tq - Tw)\ wall 


(1.5-5) 


The definition of the displacement thickness requires the mass 
rate of flow of inviscid fluid past a body thickened by an amount 6 
(measured in the N coordinate) to be equal to the viscous mass flow 
past the real body. This definition yields the following relation 
for the displacement thickness 


/ • 00 ^ f <x> 

dN = / rjpu dN , (1*5-6) 

6 ^ o 

where the displacement body is located at n = eb . When written in 
similarity variables, equation (I.5-6) becomes 

where 6 is the displacement thickness measured in the t) coordinate. 
The relation between the transformed and physical displacement thickness 
is given implicitly from equation (l.U-2) where it is seen that 6 can 
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be -written as 


Ugroj 


( 1 . 5 - 8 ) 


1. 6 Reduction to the Incompressible Case 

For the incompressible case the density is assumed constant and 
the variations in temperature are asstimed small so that the energy 
equation is iincoupled from the momenttun law. Thus, for the incompres- 
sible case p = pe = M-g = 1 in accordance with the nondimensional 


scheme and the governing equations become 


the CONTINUITY EQUATION 


^ + 2t — + F= 0, 


( 1 . 6 - 1 ) 


and the S -MOMENTUM EQUATION 


8 jR^^rdF 


^(F - 1) h^ ) 

^ J ; 


P(p2 - 1 ) + 2^F ^ + V + 231 + 2| 

dTi 

( 1 . 6 - 2 ) 


R = r/r^. 


(l.6-5a) 


h = 1 + e A dT), 

Ue^o^ ^o RJ 


(1.6- 3b) 


85 h 
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33 = 2i£He = 2L_^, 

Ug 


(i.6-5d; 


and 


I = h"- / (1 - F^) dTi . 

h5 


(l.6-3e; 


The transformed similarity coordinates | and rj are given by 


^ = ./ Ugro^^ ds , 
o 


(l. 6-ha) 


and 


Ue^oj rJ 

Ti = ~ / — dN . 

i2^ Jo h 


(l.6-4b) 


Expanding the LHS of equation (1.6-2) we obtain 


LHS 


= h2 - N + 2j (1 - F) , 

6t] \ h'^ 6 t] / h \ h h ^ / 


which can further be simplified by noting that 


(1.6-5) 


6h _ 6h 5 n _ eK V 2 I h 
Sii 6 W St) UgroJ Rd 


( 1 . 6-6a) 


Thus, hy^T^/h^ can be written as 




( 1 . 6-6b) 
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With this simplification the s-momentum equation becomes 


h2 A 


h^ ySri 


= b(F^ - 1) + 2|F ^ + V ^ + 2pl + 21 ^ 

SI OTi 0| 


(1.6-7) 

To determine the pressure in the boundary layer we must first take the 
limit of equation (l.3-8g) as the flow Mach number goes to zero. 

To do this, we write equation (l.3-8g) in terms of the inviscid 
Mach number Mg to obtain 


Pm " ^e 


1 + 


-y _ 1 2 ^ 

^ 1 1 - 




7-1 


( 1 . 6 - 8 ) 


which when expanded for small Mg before letting Mg 0 yields 


Pm = Pe I ~ Mg 0 . 


(1.6-9) 


Thus, the pressure at any point in the incompressible boundary layer 
is given by 

P(^,T1) = Pe + I Ug2/l - 

Where Uj^ is given by equation (l. 5-8e). The boundary conditions 
are given by equations (l.^-l8a) and (l.4-l8d-e). 


h2 


2~ 
+ Um I 


(1.6-10; 


■''The remaining third-order term in equation (1.6-5) has been 
neglected. 
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The incompressible forms of the boundary- layer properties are 
given as the modified skin friction coefficient 

I pSUS^ Ui V^l/wall 

where 


* 

’’’wall 


= 



wall 


(1.6-llb) 


and 

Cf = eCf . (1.6-llc) 

The displacement thickness from equation ( 1 . 5 - 7 ) becomes in the 
transformed coordinate 


6 = 



(1 - F) dri , 


and the physical coordinate value is given as 


(l.6-12a) 



dW 



R^u dN 


(l.6-12b) 


Thus, we have shown that the seconds order incompressible boundary- 
layer equations presented by Davis, Whitehead, and Wornom (197O) are 
the limit form of the equations presented in Section l.k as 0 . 
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In Chapter II, the incompressible second-order boundary-layer 
equations will be employed to study the effects of longitudinal 
curvature on wall properties for flow around a circular cylinder at 
moderate Reynolds nximbers where the displacement-body effects are 
approximately known. 



II. LONGITUDINAL CUEVATUEE AND DISPLACEMENT- 


SPEED EFFECTS ON INCOMPRESSIBLE BOUNDARY PAYEES 

2. 1 Introduction 

The major diffictilty in applying the equations formulated in 
Chapter I to a particular flow problem arises because at the present 
time, the inviscid flow parameters necessary to match the boundary- 
layer flow to the inviscid flow outside of the boundary-layer cannot 
be easily calculated to second-order for an arbitrarily shaped body. 

The purpose of this chapter is to develop a procedure for numer- 
ically solving the incompressible equations valid to second-order and 
to apply the procedure to cases where displacement- speed effects can 
be approximately determined. This application will give a test of 
curvature effects and will indirectly determine whether or not 
displacement-speed effects are Important. 

The particular flow problem to which the theory will be applied 
is that of flow around a circular cylinder. This problem was chosen 
because the effects of curvature should become important at moderate 
to low Eeynolds n\imbers and because experimental and numerical data 
in the form of wall pressure and shear stress distributions are 
generally available in the open literature. 

2. 2 Governing Equations 

The governing equations have been formulated by Davis, Whitehead, 
and Wornom ( 1970) and as was shown in Chapter I, their resulting 
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equations are the incompressible counterpart of the second-order 
compressible boundary-layer equations formulated there. 

The inviscid pressure gradient parameter p given in equation 
(l*^- 15 ) monitors the pressure gradient the inviscid flow produces in 
the I direction. It will prove beneficial later in the incompressi- 
ble study to also examine the flow in terms of what pressure gradient 
the viscous flow sees. To do this we need to look at the viscous 
pressvire gradient in the ^ direction. Differentiating equation 
(1.6-10) with respect to ^ we find that 


~ d? ® "dT 


^1. 

h 



(1-21) + 


h2 


(2.2-la) 


where 


I = h^ / (i _ F^) dq . (2.2-lb) 

^ ^ h3 


Note that 


£Pe 

dt 


+ U 


dUg 

® d^ 


= 0 . 


( 2 . 2 - 2 ) 


Rewriting equation (2.2-la) in terms of the inviscid pressure gradient 
parameter p, we find 


Sp 

65 



p(l - 21) - 2^ — 
6| 


(2.2-3) 


+ 


This wlll.be shown in Section 2.3a to be valid to second- order. 
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If we define the viscous pressure gradient parameter as 


Pv 


u| 36 


(2.2-4) 


it can be seen that the viscous pressure gradient parameter wlH 

always be opposite in sign of the actual | pressure gradient. 

Substituting equation (2.2-3) into equation (2.2-4) we obtain 
the following relation for 


“ 21 )- 2 ^ ^ ( 2 « 2 “ 5 ) 

which gives a direct measure of the longitudinal pressure gradient at 
any point in the viscous flow region. 

Wow examining the governing equations given in Section 1 . 6 , it 
can be seen that to numerically solve a general flow problem, it is 
necessary in general to know P^U), U^U), k(^), ro(|), 

cos (^), and an € . However, for the case of flow around a 

circular cylinder j = 0 and, therefore, the only necessary parameters 
are Pe(l), Ue(|), k(|), and €. 

In the present study, putting e = 0, reduces the governing equa- 
tions to the classical or first-order bo\md ary- layer equations. For 
this case Ug , the inviscid surface speed, can be taken from potential 
theory or it can be calculated from experimentally measured wall pres- 
sures using Bernoulli's equation under the assumption that the pressure 
is constant across the boundary layer. 
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However, if the governing equations are to be valid to second- 
order, then Ug must include displacement-body effects. At the pre- 
sent time this cannot be done directly. The next section seeks to 
show how this difficulty can be overcome if the wall pressure distri- 
bution is known. 

2. 3 The Pressure Gradient Parameter, Pg 

a. General Discussion .- From equation (1.6-10) we find that the 
wall pressure is given by 

Pwall = = Pe + U| I, (2.3-1) 

where 


^(l-F2)dTi. (2.3-2) 

The inviscid pressure, P^, and velocity, Ug, can be related to 
each other by first noting that to second order in the inviscid flow 
region the s component of the Navier-Stokes equation is (see 
Van Dyke (1962a)) 



u + V 8 (hu ) 6 p _ 

8s 8n 8s 


(2.3-3) 


V 


On the body surface 


Vg(s). 


U = Ue(s), 


h = 1, 


P = Pg(s), and 


Therefore, on the surface, equation (2.3-3) becomes 
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u 

e 



V 


S(hU) 


"Sir 


n=0 


dPe 

ds 


= 0 . 


(2.5-^) 


In the inviscid flow region the flow is assumed irrotational. 

This leads to the result that 

^ _ 5(hU) . (2. 3-5) 

Ss an 

Substituting this result into equation (2.3-^) we find that 





(2.3-6) 


But from Van Dyke (1962a) the inviscid normal velocity component 
can be expanded as 






+ 



On any surface without injection 


Ve = 0 + eVg^(s,o) , 


which leads to the result that 








“cTs" 


(2.3-7) 


(2.3-8) 


is third order. 
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Therefore to second- order 

(2.3-9) 


ds 


+ Ife 


ds 


= 0 , 


or 


P + i = constant . 
e 2 6 


(2.3-10) 


Evaluating the constant at the free- stream conditions gives 


P 


e 


+ 


iu2 



(2.3-11) 


Upon substituting Ug from equation (2.3-11) into equation (2.3-1) 
we obtain 


Pe 


Pwall - (1 * 2P„)I 
1-21 


(2.3-12) 


where I is given by equation (2.3-2). For classical or first-order 
boundary- layer theory e = 0 and therefore 1=0. Thus from equa- 
tion (2.3-12), one recovers the well known result that to first order, 
the inviscid surface pressure is equal to the wall pressure. To obtain 
the inviscid surface pressure valid to second order, equations 
(2.3-12), (2.3-2), and (l.2-6b) show that the actual wall pressure 
must be corrected for longitudinal curvature effects. 



Since the inviscid parameters necessary for a numerical solution 


are p and U , it' is convenient to look at equation (2.5-12) in 

t- 0 

terms of Ug. Using equation (2.3-H)> equation (2.3-12) can be 
written as 

^ •■' 1/2 

U =\ Pwall ’ , (2.3-13) 

® ; 1-21 ! 


where 


^Pwall ^^Pwall 


P ) • 


CX3 


(2.3-1^) 


Equation (2.3-13) gives the inviscid surface speed to second 
order that would he necessary to produce the observed wall pressure 
distribution. Thus, equation (2.3-13) should include the displacement 
speed. 

C_ can be computed using experimentally measured wall pres- 

■t'wall 

sure distributions, thus leaving only the integral I, (defined in 
equation (2. 3-2)) to be determined from the boundary-layer solutions. 
Assuming that the governing equations have been programmed for numeri- 
cal solutions, the integral I can be determined in the following 
manner. As was noted earlier, once the program has a Pp(|), U„(^, 

k(|), and an e , it can numerically calculate a boundary-layer solu- 
tion at a particular station. To determine I and the value of 
to second order, an iterative approach will be used. We will guess 
a value of Ug at a station, from this determine a Pg and then 
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compute a bound ary- layer profile. Taking the value of I from this 
boundary-layer profile, a new value of will be calculated from 

equation (2. 3-13) • Using this value of Ug, and new Pg will be 
determined and another iteration will be made on the boundary- layer 
solution. This process of updating the value of Ug after each 
iteration on the boundary- layer solution will be continued until I 
(and therefore Ug) converges. 

b. Solutions Aft of a Stagnation Point . - The initial guess of 
Ug at any station downstream of the stagnation point is obtained 
from the following relation^. 


U. 


= U. 


'(New Station) (Old Station) 


\ / 


As 

(New Station) 


(2.5-15) 

In this relation As is the step size in the s -coordinate and 
"ofi^New Station obtained by a linear extrapolation of the pre- 

vious two values back along the body. With the initial first guess of 
Ug and at the new station, is calculated using equation 

(l. 6-4a) and the definition. 


Since equation (2.3-15) is a two- term Taylor series expansion 
for Ug jjg^j the value of dUg/ ds used should have been the value at 

the previous station. However, for purposes of programming simplicity 
and the fact that equation ( 2 . 3 - 15 ) only supplied an initial guess, 
(dUg/ds)New station instead. 



^3 


'-’e 


(2.3-16) 


With these values of Pg, Ug, and | a bound ary- layer profile is 


computed. After each iteration the inviscid surface speed is updated 

dU. 

is updated after each iteration from the 


e 

using equation (2. 3-13)* 


ds 


following relation. 


dUe (present iteration) - Ug (previous station)^ 3 I7) 


ds 


As 


dU 

Once Ug, I, and have been updated, Pg is updated and another 

iteration is made on the solution. This process of updating the values 
dU 

of U and is continued until some criteria of convergence has 

ds 

been reached. If Ap„ < 10"5 and AF^, < 10 "^ where 

« — 'iw - 


“ Pg(last iteration) - Pg(present iteration )[, (2.3-I8) 


and 


AF 


\ 


^ /wall 


(last iteration) - (present iteration] 


(2.3-19) 

the solution was considered to have converged for this study. Note 

here that for numerical solutions of the first-order and second-order 

problems AF was held at 10 “^. 

*w 
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Solutions in a Stagnation Point Region . - At the stagnation 
point I and Ug are zero. Thus it can be seen from equation (1.6-1) 
and (1.6-7) that the | derivatives drop out of the governing equations 
which then reduce to ordinary differential equations dependent ohly 
upon the parameter p and the scale factor h. p is computed using 
equations (l.6-3c) and (2.3-16) while h is computed from equation 
(l.6-3b) using an initial guess of h to evaluate the integral for 
the first iteration and the calculated value for the next iteration and 
so forth until the boundary-layer solution converges. Because Ug and 
I are zero at the stagnation point, equation (2.3-16) which is used 
in computing p and equation (l.6-3h) which is used to compute h 
become indeterminate at the stagnation point. To obtain the values 
of Pg and '/Syue (the indeterminate term in equation (l.6-3b) at 
the stagnation point, we must expand Ug in a Taylor series about the 
stagnation point for small s. Thus Ug(0 + s) becomes 

Ue = u;s + ... , (2.3-20) 

where U^ is the value of dUg/ds at the stagnation point. 

Using equations (2.3-20), (l.6-4a),and (2.3-16) we find that at 
the stagnation point. 


Pe = 


(2.3-21) 


and 



/ 
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( 2 . 3 - 22 ) 


Thus the major difficulty at the stagnation point will be determining 
the proper value of . To do this, note that from equation (2-3-13) 
we see that for Ug to be equal to zero at the stagnation point, 
must be 1 there. Both Homann (1936) and Grove, Shair, Petersen, 
and Acrivos (1964) have shown this to be true to second order. 

Expanding C about the stagnation point, we obtain 

-t’wall 


^wall 


1 + AoS”^ 


(2.3-23) 


Substituting the expansions of Ug and C about the stagnation 

^wall 

point into equation (2-3-13) and taking the limit as s goes to zero, 
we find that 


"^Homann (1936) and Grove, Shair, Petersen, and Acrivos (1964) have 
shown analytically from a boundary-layer analysis that at the stagna- 
tion point 


C = 1 

^wall 



(2.3-22a) 


where a is a constant. Using a boundary-layer concept, they found 
a to be equal to 8. More recently, Takarai and Keller (1969) assumed 
equation (2. 3-22a) to be correct and estimated a to be 5.985 through 
comparison with their numerical solutions of the Navier-Stokes equations 
at Reynolds numbers of 50 and 60. Dennis and Chang's (1969) numerical 
solutions for Reynolds numbers from 50 to 100 also give a ~ 6. 


k6 


u; = iAaJ'u - 


(2.3-2h) 


where is the value of the integral I at the stagnation point. 

Noting that for a circular cylinder k(|) =1, equation (2.5-2) 
in boundary-layer coordinates becomes 


I 


o 



(2.3-25) 


As an initial guess of at the stagnation point, a second-order 

estimate of equation (2. 3-2k) will be made. To do this, we only need 
a first-order estimate of the integral in equation (2.3-25) and this 
can be found by substituting the following expansions for u and Ug, 

u = up + eu2 + ... (2.5-26) 


= U, 


ei 


+ eUc 


(2.5-27) 


into equation (2.5-25) and keeping only terras of first order in that 
integral. Therefore equation (2.5-25) becomes 


I 


o 




(2.5-28) 


which can be expanded into 



4? 



(2.5-29) 


¥e can now identify the two integrals in equation (2.3-29) as 
"being the first-order displacement thickness 6p , and the first-order 
momentum thickness so that we may write that 


- e(S-i + 0n ) . 


(2.3-30) 


From Schlichting (1968), the first-order values of Sp ®1 

a stagnation point for a cylinder are 

&^=0.ij-58, (2.5-31) 

and 

0p = 0.206, (2.3-32) 

where the value of dU^/ds at the stagnation point used to obtain 
these values was taken from potential theory. 

Thus the initial guess of at the stagnation point will be 

taken from the following relation 

Ui = (1 - 1.328e)-l/2, (2.3-35) 
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where = 0.664e has been substituted into equation (2.3-24). The 
constant is to be determined from the wall pressure distribution 

in accord with equation (2.3-23)- 

After each iteration on the boundary-layer solution, is 

updated using equation (2-3-24) until s.atisfies the following 

convergence criteria Wq < 10“ ^ where 

ZUq = |uQ(present iteration) - UQ(last iteration)] , (2-3-34) 

and Zli' < 10“ 5 • was defined in equation (2.3-19)- 

'Iw — % 

Table I gives the initial values of Uq obtained from equation 
(2-3-33) using a series of various wall pressure distributions that 
will be discussed in detail in Section 2.5b. Also included in this 
table are the final values of obtained from the iterated boundary- 

layer solutions. 

2.4 Outline of the Nunerical Solution Technique 

The second-order boundary- layer equations have been programmed 
for numerical solution using an implicit finite-difference method 
similar to that developed by Fliigge-Lotz and Blottner (1962) and as 
applied to higher-order boundary- layer theory by Davis, Whitehead, 
and Wornom (1970).'*' We want to outline here the procedure used to 
numerically solve the governing equations when correcting the wall 
pressure for curvature effects to obtain the second-order inviscid 

"''Details of the method of solution are given in Appendix A. 
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surface pressure. For a general station away from the stagnation 
point the following procedure is en^loyed: 

1. The initial estimate of Ug is taken from equation (2. 3-15) • 
The initial estimate of dU^/ds is made by extrapolating the previous 
two station values to the new station. 

2. I and Pg are computed using equations (l.6-4a) and (2.3-16) 
and one iteration on the boundary-layer solution is performed. 

3" Using equations (2. 3-18-19) ^ check is made to see whether 
the solution has converged. If not, the value of I from this solu- 
tion is used to update the values of Ug and dU^/ds through equations 

(2.5-13) and (2.3-17)- 

U. Steps 2 and 3 are repeated until the solution converges. 

At the stagnation point the initial guess of dU^/ds is taken 
from equation (2.3-33)- For each iteration thereafter, dU^/ds was 
updated using equation (2.3-24). 

2. 5 Discussion of Results 

a. Wall Pressure Data. -Four sets of experimental wall pressure 
data were used in this study. They were the experimental data of 
Thom (1953) at Re^ = 17^, and three sets of experimental data by 
Grove, Shair, Petersen, and Acrivos (1964) -- two at Re^j = I75 and 
the other at Re^j = 177- The Reynolds number Re^ is. defined as 
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+ 


Red 


U* d-^ 




* 

00 




(2.5-1) 


where d* = diameter of the test cylinder. 

The d*/h* ratios employed in these experiments were 0.025^ 0.1, 
0.05, and 0.05 respectively, where 


d* _ diameter of test cylinder 
h* width of test section 


All four sets of experimental data were assumed to represent 
steady flow conditions althoiigh only the data at Reynolds numbers I 75 
and 177 were taken using a splitter plate behind the cylinder to insure 
wake stability. Figiire 2.2 shows the variation of the experimentally 
measured wall pressure distributions for Reynolds numbers I74, I75 
(d*/h* = 0.05 distribution only), and I77. The exact reasons for the 
wide variations is not known but it seems to closely associated with 
the splitter plate influence. 

Two numerically obtained wall pressure distributions were also 
used in this study. The first was obtained by Kawaguti (1955) from 

a numerical solution of the steady-state incompressible Wavier-Stokes 

"^The second-order bound ary- layer equations were nondimensionalized 
with respect to the nose radius of curvature shown in Figure 2.1. 
Therefore, the Reynolds numbers appearing in the governing equations 
through e is based on the radius of the cylinder. However, when 
referring to the different Reynolds number cases, the Reynolds number 
stated will be based on the cylinder diameter which is the Reynolds 
number used in reporting the experimental and numerical data. 
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equations at Reynolds number 4o. The second one was obtained by Thoman 
and Szewczyk ( 1969 ) from a numerical solution of the time-dependent 
incompressible Navier-stokes equations at Reynolds number 200. 

b. Curvefit Representation of the Pressure Distributions . - In 
order to be able to compute the surface pressure at arbitrary stations 
along the surface, least- square curvefits were applied to all the wall 
pressure distributions used here. 

Since steady- flow conditions are asstimed, the wall pressure 
should be symmetrical with respect to the body coordinate measured 
from the stagnation point. For this reason, the four sets of experi- 
mental data and the numerical wall pressure distributions were fitted 
to a curve over the interval -l 80 ° to +l 80 ° by a least- squares technique 
using the ratio of two Chebyshev polynomials of the form 


^wall 


^wall ~ 

1 p*U*^ 

2 00 00 


IP 

V C T (x) 
^ n n 


) ^n'^n ^ ^ ) 

n=t) 


(2.5-3) 


where IP is the number of terras in the numerator expansion and IQ 
is the number of terms in the denominator. given by the 

recursion relation 


Tn(x) = 2xTn_i(x) - 1 ^. 2 (x) , 


( 2 . 5 -^) 
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where Tq - 1, = x, T2 = 2x^ - 1, . . . and x is contained on the. 

interval -1 < x < 1. Further details concerning Chebyshev polynomials 
are given by Carnahan, Luther, and Wilkes (1969)- Table II lists the 
coefficients C^, C2, Ci|,, . . . and D^, D2, D^, . . . for the different 

curvefits used. Because the wall pressure was taken to be an even 
function of the body coordinate s, the odd coefficients were deter- 
mined to be zero when the data was fitted over the interval -l8o° to +l80°* 

The criteria upon which the selection of the curvefits used in 
this study was based will be presented in Section 2. 5h where some of 
the curvefit difficulties encountered will also be discussed. 

The wall pressure data was not tabulated in the literature for 
the Reynolds numbers ^0, 175^ 177^ ^•i^d 200 cases. Table III gives the 
data points for these cases as they were read from graphs given in the 
cited references. These points were the basis for the wall pressure 
distributions employed in this study. Figure 2. 3a~f show comparisons 
between the pressure distribution curvefits selected for this study and 
the data points used in obtaining such. 

It was pointed out in Section 2. 5c that at the stagnation point 
^Pwall equal to 1 to second order. Note that this is not 

true of the wall pressure data shown in Figures 2. 5a-d. This differ- 
ence at the stagnation point is an indication of third and higher-order 
effects present in the data. 

Figure 2. 3e and 2. 3f show ^ to be approximately 1 at the 

stagnation point for the Reynolds number 200 wall pressure data. This 
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appears to be an inconsistency as equation (2. 3-22a) implies that the 
stagnation value should be approximately 1.05 or 1.04. It is not 
known how serious this inconsistency is, but some difficulties were 
encountered in trying to find accurate curvefit representations of this 
particular wall pressure data in the stagnation region (see Sect. 2.5h). 
This inconsistency may have been the cause, but this is not known for 
sure. 


At other points along the body away from the stagnation point, the 
third-order corrections influence on the wall pressure distribution is 


not known. For want of a better method, the stagnation point values of 


these higher-order effects (i.e., a/Re^) were assumed representative 


for all s and subsequently subtracted out of all the pressure data. 
This amounts to translating the wall pressure curvefit distribution 


after the data has been curvefitted"*" so that 
stagnation point 


Pwall 


equals 1 at the 


When the wall pressure data was translated before ciirvefitting, 
the resulting curvefit still required shifting because it usually was 
not identically equal to 1 at the stagnation point. Such an unequal, 
shifting of the data and the curvefit resulted in a poorer representa- 
tion of the wall pressure data than the above method. 

’’"‘'Other authors, namely, Diraopoulos and Hanratty (1968) fitted the 
Reynolds number 174 wall pressure data to a polynomial of the form 


p* 

■^wall 


p* 

00 


r 

2 




(2.5-5) 


^ Using this pressure distribution, they obtained a Blasius series solu- 
tion of the classical boundary-layer equations for the shear-stress 
distribution. It should be noted that in doing their Blasius solution, 
Dimopoulos and Hanratty (1968) did not need to shift their pressure 
distribution since the stagnation value of does not appear in 

the final Blasius shear- stress equation. 


c. Effects of Curvature and Displacement Speed on Wall Pressure 
Distributions . - Figures 2. 4a-e show the following wall pressure distri- 
butions for all of the Reynolds number cases studied. Curve (l) gives 
the inviscid surface pressure distribution obtained from potential 
theory and it is the same in all the figures. Curve (2) shows the 
wall pressure distribution taken from either experimental or numerical 
data. Curve (3) shows the inviscid surface pressure valid to second 
order obtained from equation (2.3-12) and the second-order boundary- 
layer solution. Looking at the results in a different light, curve 
(3) shows the inviscid surface pressure to second order that would be 
required to produce the observed wall pressure distribution shown by 
curve (2). As such, the difference between curve (3) and curve (2) 
can be attributed to curvature through equation (2.3-12). Without the 
curvature correction, eauation (2. 3-12) implies the wall pressure is 
equal to the inviscid siirface pressure. 

Fig\ires 2.ha-e also give some indication of the effect the 
boundary-layer displacement thickness has on the inviscid surface 
pressure. If one assumes the inviscid surface pressure distribution 

"^Due to the presence of a wake behind a cylinder even at high 
Reynolds numbers, the inviscid potential pressure distribution is 
meaningless for 0 > 100°. 

’*""''Curve(4), shown in Figure 2.4c only, shows the wall pressure 
that results from a solution of the second-order equations using the 
potential theory inviscid surface pressure distribution to compute 3g. 
Thus, curve (4) gives the wall pressure that would exist if the 
inviscid pressiire were given by curve (l). 
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at infinite Reynolds number to be given by potential theory, ^ then the 
difference between cui’ves (l) and (3) of Figures 2.4a-e can be viewed 
as a measure of the second- order contribution to the inviscid flow 
(i.e., that due to flow over the displacement body). 

d. Displacement- Body Effects . -For the governing equations to be 
valid to second order, the inviscid surface speed Ug must be the 
inviscid surface speed on the cylinder surface due to flow over the 
displacement body. In Section 1.4 we defined a displacement thickness 
S such that the mass rate of flow of inviscid fluid past a body 
thickened by an amount 8 is equal to the amount of viscous flow past 
the real body. This definition yielded the following relation for the 
displacement thickness 5 



In terms of the displacement thickness 5, the displacement body for 
a cylinder is then given by 1 + eS. 

Figures 2.5a-e show the displacement bodies obtained in this study 
for solutions of the complete"*^ second-order equations. It should be 
noted also that for the complete second-order solution, the inviscid 
surface pressure shown in Figures 2.4a-e is the inviscid surface 

"^Due to the presence of a wake behind a cylinder even at very 
high Reynolds numbers, this point is uncertain. 

++ 

Complete here refers to solutions of the second-order equations 
where the inviscid flow parameters are also valid to second-order. 
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pressure which boundary- layer theory says is induced on the cylinder 
surface due to flow over these displacement bodies. The thickness of 
the boundary layer for the Reynolds number 4o case gives an indication 
that these solutions are an extreme test of the limitations of the 
boundary-layer concept. 

Figure 2.6 shows the displacement body obtained from a classical 
boundary- layer solution at Reynolds number 177 using the inviscid 
surface pressure taken from potential theory to compute Pg. Whereas 
in Figures 2.5a-e the inviscid flow sees the approximate displacement 
body, this is not true for the solution shown in Figure 2.6 where the 
inviscid solution is also for a cylinder. Note that Figure 2.6 gives 
evidence of a singularity in the displacement thickness derivative at 
separation. However, Figure 2. 5d shows that when the viscous flow 
was made consistent with the inviscid flow (i.e., the complete second- 
order solution), there appeared to be no evidence of such an impending 
singularity at separation. 

e. The Viscous Pressure Gradient Parameter Py. - Figure 2.7 
shows the values of the inviscid and viscous pressure gradient param- 
eters on the cylinder surface for the Reynolds number 177 case."*"^ 

Also shown are the values at the outer reaches of the viscous flow 
region (i.e., where F = 1.0000). The large variation in the viscous 

’’’Brown and Stewartson (1969) have shown that if there is a 
singularity in the shear stress at separation, there will also be a 
singularity in the derivative of the displacement thickness. 

^’’’The other Reynolds numbers cases were similar. 



pressure gradient parameter is an indication of large pressure gradi- 
ents existing across the boundary layer. From Figure 2. 7 and also 
equations (l. 4 -l 8 e), (2.2-lb), and (2.2-5) it can be see that 

Pv ^ 3 as N CO . (2.5-7) 

It is often suggested that separation of the classical boundary-layer 

flow from a body results from the boundary layer being subjected to 

a positive pressure gradient. For classical bound ary- layer theory one 

need only monitor the inviscid pressure gradient parameter Pg to 

determine when the boundary layer sees a positive pressure gradient. 

The reason for this is that for classical bound ary- layer theory the 

pressure is assumed constant across the boundary layer and as a result, 

the viscous pressure gradient in the | direction is equal to the 

Inviscid pressure gradient. However, for the second-order equations 

the viscous pressure gradient is not equal to the inviscid pressure 

gradient. For this case the viscous pressure gradient parameter p^ 

+ 

must be monitored. Monitoring p^. Figure 2 . 7 shows the portion of 
the boundary layer nearest the wall undergoing a positive viscous 
pressure gradient at separation. This causes the flow to separate at 
the wall and in turn causes the rest of the boundary layer to separate 
with it. 


Note from Figure 2.7 that if one were to monitor the inviscid 
pressure gradient parameter p, he would led to the conclusion that the 
boundary layer was undergoing a negative pressure gradient at separation. 
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f- Shear Stress Distributions .- Figures 2. 8a-e show the shear 
stress distributions obtained in this study. In each of these, 
curve (l) shows the solution to the equations valid to second order 
using the inviscid surface conditions obtained from equation (2.3-13) 
by iteration and as such, represents the complete second-order solu- 
tion. The remaining curves will now be discussed in terms of how they 
compare with curve (l). Comparison of the complete second-order 
solution, curve (l), with experimental results will be left for a 
later section. 

Curve (2) shows the solution to the classical first-order boundary- 
layer equations using the inviscid surface pressure distribution 
obtained from potential theory. This solution is the same for all 
Reynolds numbers. The shear stress given by this solution can be 
seen to be very much higher than curve (l). 

Curve (3) shows the solution to the classical bound ary- layer 
equations taking the inviscid surface pressure to be approximately 
equal to the observed wall pressure distribution. It was quite 
surprising to find this solution virtually equal to the second-order 
solution except near separation. This result can be rationalized 
in the following way. Wear the wall where N is small, h is 
approximately equal to 1, and u is approximately equal to zero. 
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Substituting these approximations into the governing equations gives 
a zero normal pressure gradient and reduces the remaining two equations 
to the classical bound ary- layer equations. Therefore, the flow near 
the wall is approximately governed by the classical boundary-layer 
equations even when curvature terms are included. Since both of the 
solutions, curve (l) and curve (3), see the same wall pressure gradient, 
this apparently results in their shear stress distributions being 
approximately equal. We note in passing that, based on this point, 
it is natural to conclude that an accurate means of obtaining the 
wall pressures is needed in order to calculate the wall shear. 

Curve (4) (shown in Figure 2.8c only) gives the solution to the 
second-order ^equations using the inviscid surface pressure distribu- 
tions from potential theory to compute (3^. This solution makes no 
attempt to correct the inviscid surface pressixre (potential) for 
curvature effects. The wall pressure for this solution is not the 
observed wall pressure distribution but the wall pressure that comes 
out of equation (2.3-1)- For the Reynolds number I 75 case, this wall 
pressure is shown by curve (4) in Figure 2.4c -- the wide difference 


"''The second-order incompressible boundary- layer equations are 
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between this and the actual wall pressure apparently causing the 
difference in the shear stress distributions observed in Figure 2.8c. 
Note that for this case the boundary- layer solution failed to con- 
verge after 30 iterations "*"( the maximum allowed for this particular 
case) to the degree of accuracy required for this study. It is also 
worthy of note that this solution lies closer to the complete second- 
order solution than curve (2) shown in Figure 2.8c. 

Curve (5) shows the solution to the second-order equations using 
the wall pressure distribution to conq)ute ,8g, that is^ taking the 
inviscid surface pressure to be approximately equal to the wall 
pressure distribution. This assumption is often made since the measured 
wall presstire distribution should contain displacement effects. If the 
pressure is approximately constant across the boundary layer, then it 
would seem that taking the observed wall pressure distribution to be 
approximately equal to the inviscid surface pressure would be a valid 
assumption. However, as was noted earlier, this assumption is not 
a very good one for the Reynolds numbers studied here as can be seen 
from the results in Figures 2.4a-e and also from the shear stress 
results shown in Figures 2. 8a- e. 

g* Comparisons with Experimental and Numerical Shear Stress Data . - 
In order to properly compare the numerical shear stress distributions 
with experimental data, wall pressure and shear stress distributions 

"'’Thirty iterations at that particular station. 

■‘"'■CorTrputer output showed that the solution was in the process 
of converging when the 50 iteration limit was reached. 
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were needed which were taken at the same Reynolds number and under 
the same test conditions. However, none were to be found in the 
literature. 

The experimentally measured wall pressure and shear stress dis- 
tributions appear to be very sensitive to the test parameters d*/h*, 
defined in equation (2.5-2) and _c*/d* defined as the distance from 
the cylinder centerline to the front edge of the wake splitter plate 
ratioed to the diameter of the cylinder. 

The influence of c*/d* ratio on the measured wall pressure can 
be partly seen in Figure 2.2. For a Reynolds number difference of 2 
between the 175 and 177 case, there is a slight change in the wall 
pressure data. For the 175 case, c*/<i* was approximately 2.5^ 
d*/h* was 0.05. For the Reynolds number 177 case, c*/d* was 5* ^ 
and d*/h* was 0.05- 

The effect of the d^/h* ratio on the measured wall pressure 
distribution was shown experimentally by Grove, Shair, Petersen, 
and Acrivos (1964). To illustrate here how much the d*/h* ratios 
influence the numerical shear stress results, the measured wall 
pressure distribution at Reynolds niimber 175 with c*/d* = 2.5 and 
d*/h* = 0.1 was employed in addition to the Reynolds number I75 case 
previously mentioned with a d*/ h* = 0.05 and a £*/d* = 2.5. Figure 
2. 9 shows a significant difference in the two Reynolds number I75 
wall pressure distributions. Figure 2.10, which was obtained from a 
solution of the classical bound ary- layer equations taking the inviscid 
surface pressure to be approximately equal to the observed wall 
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distribution"'" shows a comparison between the two resulting shear stress 
distributions. Note from Figure 2.10, that the d*/h* ratio does not 
appear to have much influence on the shear stress distribution in the 
stagnation region. However, at points further back along the body, 
the influence becomes noticeable. Figure 2.*10 strongly suggests that 
when comparing numerical shear stress distributions obtained by this 
method with experimental shear stress data that, in addition to matching 
the Reynolds number, one should also match as closely as possible the 
d*/h* ratio at which the wall pressure and shear stress data were 
taken. 

Of all the experimental shear stress data found in the literature, 
the closest (in terms of the Reynolds number and test conditions) to 
the Reynolds number I 75 case where a splitter plate was used, was the 
data at Reynolds number l 6 ? with a c*/d* ratio of 2.5 and a d*/h* 
ratio of O.O 853 . Figure 2.11 shows a comparison between the classical 
boundary- layer shear stress solution using the measured wall pressure 
to compute and the experimental shear stress data at Reynolds 

number I 67 . Since Figures 2-9 and 2.10 show the wall pressure and 
the resulting shear stress distribution to be very sensitive to the 
d*/h* ratio, this fact should be kept in mind when studying Figure 
2 . 11 . 


In Section 2. 5 f, this solution was shown to be virtually the 
same as the complete second-order solution. Therefore, because this 
solution requires less computer time to run, it will be employed for 
comparison purposes. 
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For the Reynolds number 17^ wall pressure distribution which was 
taken without the aid of a splitter plate, the nearest experimental 
shear stress data was at Reynolds number 151- However, the d*/h* 
ratio for the Reynolds number 17^ case was 0.025, whereas for the 
Reynolds number I5I case, d*/h* was 0.0835* Because of the wide 
difference between the d*/h* ratios used in obtaining the experimen- 
tal wall pressure and shear stress distributions for this case, it was 
felt that any conclusions arrived at by comparison with this data 
would seem open to question. Therefore this comparison was not made. 

Similarly, for the Reynolds number 177 case, the large differences 
in d*/h* ratios used to obtain the wall pressure and shear stress 
data prevented a comparison of such. 

Note, however, that DimopoilLos and Hanratty (1968) did compare 
the Reynolds number I5I shear stress data with a Blasius series solu- 
tion using a curvefit of the Reynolds number wall pressure by 
Thom (1933) to compute 3^. The good agreement between experiment 
and theory which they report, apparently results from an inconsistency 
in their Blasius solution. Using the Blasius series taken from 
Rosenhead (1963),"^^ and the wall pressure representations and constants 


‘‘‘Equation (2.5-5) gives the form of their curvefit representation. 

^^here seems to be some question as to which Blasius series 
solution Dimopoulos and Hanratty (l968)used. Dimopoulos and Hanratty 
(1968) report they used the Blasius solution given by Rosenhead (1963). 
However, Dimopoulos (1968) reported that he used the Blasius solution 
given by Schlichting (i960). 
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of Dimopoulos and Hanratty ( 1968 ) (see eq. (2-5-5) and Table IV) we 
coiild not reproduce their reported shear stress distributions. Figure 
2.12 shows the Blasius shear stress distribution obtained in this study 
for the Reynolds number 174 case using the curvefit coefficients of 
Dimopoulos and Hanratty ( 1968 ). Also shown in Figure 2.12 is the Blasius 
solution reported by Dimopoulos and Hanratty ( 1968 ) along with the 
Reynolds number I 5 I shear stress data with which they compared. 

Two checks were made on our Blasius solution. Taking the inviscid 
surface speed on a cylinder (U^ = 2 sin 0 ), substitution was made into 
Bernoulli's equation and expanded to obtain a pressure relation of the 
form of equation ( 2 . 5 - 5 ). Picking off the coefficients in equation(2. 5-5) 
substitution was then made into our Blasius equation which was written 
in terms of these coefficients. Our Blasius solution, using the con- 
stants taken from potential theory, compared exactly with that reported 
by Schlichting ( 1968 ). The second check was a finite- difference solu- 
tion to the same equation using the wall pressure curvefit to compute 
Pg. Figure 2.12 shows that our Blasius solution agrees with the finite- 
difference solution except near separation. 

Figure 2.12 shows that their good agreement was apparently due to 
an inconsistency in their Blasius solution. Appendix B gives the 
Blasius solution obtained in this study in terms of the constants in 
the curvefit equation given by Dimopoulos and Hanratty (I 968 ). 

Because the Reynolds numbers 4o and 200 wall pressure distribu- 
tions were taken from niimerical solutions of the full Navier-Stokes 
equations, the shear stress distributions were available from the same 


65 


"test" conditions. Using the following relation, the shear stress 
distributions were obtained from the vorticity distributions on the 
cylinder surface reported by Kawaguti (1953) and Thoman and Szewczyk 

(1969) 


- -2U , , 

Cf=— (2.5-11) 

\/rT 

Figures 2.13 and 2.14 show comparisons between the boundary-layer 
shear stress distributions and the shear stress distributions obtained 
by Kawaguti (1953) and Thoman and Szewczyk ( 1969 ). Figures 2.13 and 
2.14 shows good agreement between the boundary- layer solution and the 
full Navier-Stokes solution. Why there seems to be better agreement 
between the boundary- layer solution and the Navier-Stokes solution 
at Reynolds number 40 rather than Reynolds niomber 200 is not known. 

h. Curvefitting Diffictilties . - The selection of the particiilar 
wall pressure curvefits used in this study was based on two criteria. 
The first criteria was that when a particular wall pressure curvefit 
was used to compute 3^, the Pg distribution computed had to be a 
well-behaved function. The second criteria consisted of selecting 
(from among all the curvefits that satisfied the first criteria), the 
curvefit which to the eye, appeared to best represent the wall pressure 
data. 

Since most of the curvefit difficulties encountered were with 
the Reynolds n\unber 200 wall pressure data, this case will be used for 
illustration. Two curvefits of the Reynolds number 200 wall pressure 
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data will be discussed -- curvefit A and curvefit B. Curvefit A, 
shown in Figure 2.3e, appeared to give the best representation of the 
Reynolds number 200 pressure data of all the different curvefits 
tried. However, it did not satisfy the first criteria. Figure 2.15 
shows the Pg distribution that results when computed using this wall 
pressure distribution. The unexpected flat portion from approximately 
20° to t-0° is what is termed here as a not well-behaved function. How 
this flat part of the distribution affects the boimdary-layer 

shear stress results can be seen in Figure 2.16. Up to the flat part 
of the Pg distribution, the boundary-layer shear stress agrees fairly 
well with the more exact Navier-Stokes solutions of Thoman and Szewcyzk 
(1969). Then the effect of the flat portion of the p^ distribution 
seems to be to shift the n\omerical results away from the Navier-Stokes 
solution. The reason for the horizontal portion in Figure 2.I5 can 
be traced back to the derivative of the wall pressure shown in Figure 
2.17. Whenever the wall pressure derivative became flat in the approxi- 
mate area shown, the p^ distribution computed from that cirrvefit 
distribution also changed in the same area. Figure 2. 3s shows curvefit 
B which was selected as the best curvefit satisfying both of the cri- 
teria given above. Figures 2.18 and 2.19 show that both the pressure 
derivative and the computed Pg using curvefit B are well behaved. 
Figure 2.20 shows a comparison between the classical boundary-layer 
shear stress and the numerical shear stress. 

It should be noted here that, although the wall pressure curvefits 
selected were the best we could obtain, it was felt that some. 
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particularly curvefit B of the Reynolds number 200 wall pressure data, 
did not represent the data accurately enough. Examining Figure 2. 3e, 
it can be seen, for example, that curvefit B does not represent the 
stagnation region wall pressure data very well. Upon examining the 
shear stress results obtained by using curvefit B to coir^jute pg and 
the shear stress results obtained by using curvefit A to compute 
it is noted that curvefit A, which gave the best representation of the 
stagnation region wall pressure data also compared more favorably with 
the experimental results in that region. Further evidence that the 
numerical shear stress is sensitive to the wall pressure representa- 
tion can be seen from Figures 2.21 and 2.22. Figure 2.21 shows a 
comparison between two different wall pressure representations of the 
Reynolds number 17^ data. Both appear to be fairly good representa- 
tions of the wall pressure data. However, Figure 2.22 shows that when 
both pressure distributions are used to conipute Pg and a classical 
boundary-layer solution, one pressure distribution leads to separation 
whereas the other does not. 

2. 6 Conclusions 

In summary, the accuracy of the method presented here for studying 
the effects of curvature on boundary layers at moderate Reynolds numbers 
seems to be limited by the accuracy of the measured wall pressure data 
and the accuracy of the wall pressure curvefit representation. Although 

'*'ln Section 2.5b it was noted that the stagnation value of C 
appeared to be low (implied by eq. 2. 5-22a)). Whether this ^wall 

contributed to the poor curvefit representation in the stagnation 
region is not known. 
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there was a lack of experimental shear stress data taken under the same 
test conditions as the wall pressure data, what data there was strongly 
indicates that given a very accurate wall pressure representation, the 
agreement between the boundary-layer shear stress and the experimental 
shear stress data would be very good, even for Reynolds numbers as low 
as ^0 for cylinders. In addition, we conclude the following; 

1. For all the Reynolds number cases studied, the effects of 
longitudinal curvature on the boundary-layer wall pressure are signifi- 
cant as attested to by Figures 2. 4a-e. This implies that if one knew 
the second-order inviscid surface conditions and wanted to determine 
the wall press\ire (viscous), the omission of longitudinal curvature 
effects would result in large errors in the flow properties at the 
wall. 

2. For the Reynolds number cases studied, large pressure gradients 
(implied above) were found to exist across the boundary layer. This 
would seem to exclude any hope of obtaining a good representation of 
the flow from classical boundary-layer theory. However, to the con- 
trary, it was found that if one had a measured wall pressure distribu- 
tion and was only interested in shear stress results, a classical 
boundary-layer solution using the measured wall pressure distribution 

to compute Pg would be sufficient. 

3. For the complete second-order solution, there appeared to be 
no evidence of an oncoming singularity in the displacement thickness 
on a cylinder at separation such as is indicated by a classical 
boundary-layer solution. 



III. A STUDY OF LONGITUDINAL CURVATURE EFFECTS OW COMPRESSIBLE 
BOUKDARY- LAYER FLOW WITH DISPLACEMENT THICKKESS INTERACTION 

5 . 1 Introduction 

In Chapter II, the effects of longitudinal curvature on incompres- 
sible boundary-layer flow were studied for the case of flow around a 
circular cylinder by taking the measured wall pressure distribution 
and determining what the inviscid properties were that produce that 
wall pressure. By this indirect method the effects of longitudinal 
curvature and displacement speed were ascertained to second order. 

In this chapter, we want to study the effects of curvature on 
cofl^jressible boundary-layer flow for a case where the displacement- 
body effects can be computed directly using some approximate method to 
compute the inviscid flow properties. Attention then is drawn toward 
supersonic flows where a raiLLtitude of approximate schemes are available 
- for example, linear thin airfoil theory or other comparable models. 

For this reason, it was decided to extend the finite-difference 
work of Flugge-Lotz and Blottner (1962) in which they computed 
directly the displacement- body effects on a flat plate in h3^ersonic 
flow. Here interest will be directed toward curved body surfaces at 
hypersonic speeds. 

Although the analysis by Cheng, Hall, Golian, and Hertzberg (196I) 
was not restricted to flat plates, it was not until I969 that the 
curved surface case was studied. At that time, Stollery (1969) 
extended the analysis of Cheng et al. (196I) and Sullivan (1969) to 
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cover curved body surfaces in regions of strong and weak interaction. 
These approximate theories did not, however, take into account varia- 
tions in pressure across the displacement body or boundary layer due 
to the curvature of the body as will be done in the present study. 

In addition to studying curvature effects, this chapter also seeks 
to shed light on the possible existence of a saddle-point type behavior 
that can emerge in these types of problems from numerical solutions. 
Such behavior for the strong interaction problem was observed by 
Garvine ( 1968 ) with his viscous shock-layer flat-plate solution. 

Baum ( 1968 ) observed similar findings with his supersonic boundary- 
layer solutions. However, the more recent papers of Sullivan (I 969 ), 
Stollery (1970); and Cheng, Chen, Mobley, and Huber (1970) did not 
encounter a saddle-point behavior. If indeed a saddle-point behavior 
does exist for this type of interaction problem, it should show up 
in the present study and will have to be overcome. Once the inter- 
action process is modeled correctly, the way will be clear for further 
studies on curved surfaces in supersonic flows. 

The particular problem to be investigated here is that of two- 
dimensional laminar flow up a cubic compression ramp at Mach numbers 
6 , 8 , and 12.25- This problem was chosen because the displacement- 
body effects should be sizable and the curvature effects can be easily 
handled since the surface curvature is everywhere analytic. In 
addition there exists some experimental data due to Stollery ( 197 O) 
for the Mach number 12.25 case which will serve as a basis of com- 
parison in this study. 
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3- 2 Governing Equations 

a. Viscous Region . - The governing equations for the viscid region 
are the second- order compressible bound ary -layer equations presented 

in Section 1.4 of Chapter I. Details of the numerical method of 
solution used here are given in Appendix A. 

b. Inviscid Region . - Examining the governing equations for the 
viscous region, it is seen that the inviscid properties necessary to 
compute the boundary-layer solution at a particular station are Pg , 

Tg, Pg, Pg, Ug, dUg/d|, and i . In this section, it will be 
shown how to relate these to the inviscid pressure Pg and pressure 
gradient dPg/ds at the desired station. Once Pg, Pg, and Ug are 
known, | is determined by mamerical integration employing a trapezoidal 
r\ile. 

For this study, the tangent-wedge formiila is employed to compute 
the inviscid pressure on the displacement body under the assimption 
that the displacement body is a streamline for the inviscid flow model. 
With this approximation the pressure at any point on the displacement 
body is taken to be equal to the pressure on a wedge whose half -angle 
equals the local inclination angle of the displacement- body tangent. 
Physically this approximation gains its validity from the fact that 
at hypersonic speeds the shock layer is sufficiently thin on slender 
bodies to assume that the variations in pressure and streamline angles 
will be small across the shock layer. Thus, the surface values are 
approximately the same as those at the shock. In addition to the 
above limitations, the tangent-wedge approximation can be seen to 
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neglect centrifugal- force effects that result from curvature of the 
body. None the less, this pressure law has become generally acceptable 
as a model of hypersonic inviscid flows and will be employed here due 
to its extreme simplicity. It was also adopted by Fliigge-Lotz and 
Blottner (1962), Sullivan (1969)^ Stollery (1970 and other authors 
and thus its use here allows a direct con5>arison with other viscous 
solutions. Another important feature as pointed out by Sullivan (1969) 
is that the tangent-wedge approximation gives a single relation valid 
for both the strong and weak interaction regimes. 

From Hayes and Probstein (1959) the tangent-wedge pressure law 
can be written as 
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(5.2-la) 


where 

9 = ©wall + eg, (5.2-lb) 

and it is understood that the local inclination angle, P, of the 
displacement-body tangent is sufficiently small as to make the 
sin 0 ~ 0 a valid approximation. 

In equation (5.2-la), is the inviscid pressure induced on 

the displacement- body surface. To determine the inviscid pressure on 
the actual body in terms of the pressure on the displacement body, we 
first expand the pressure about the body surface in a Maclaurin series 
thus obtaining 
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P(s,o + eS) = 



( 3 . 2 - 2 ) 


Noting that P(s,o) = Pg, equation (3.2-2) can be solved for Pg 
to obtain 


P 


e 



(5.2-3) 


From equation ( 3 . 2-3)j it is seen that to obtain Pg formally valid to 


second order we need only to retain terras in the j terra to 

s, o 

0 ( 1 ). To do this, we note that the inviscid n-moraentum equation 
written in terms of the present coordinate system as given by Van Dyke 

( 1962 c) is 


» = -p 2 ? v„ . ink 

Bn \h Bs h 


( 5 . 2 - 1 ^) 


If we now apply the ordering scheme given by Van Dyke ( 1962 c), 
that is. 


U = U3_ + eUp + . . . 


p = pp + €p 2 + ... 


V = Vp + eVp + . . . 


(5.2-5a) 


(5.2-5b) 


(3.2-5c) 


to equation (5.2-4), we obtain 
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(3.2-6) 


On the body surface h = 1 and V2_(s,o) = 0, thus to 0(l) equation 
(3.2-6) becomes 


|, = Kp 3 _(s,o)uf(s,o) . (3.2-7a) 

s, o 

Noting that pj_(s,o) and U]_(s,o) are sin^ily the inviscid surface 
values for the case where displacement- body effects are neglected^ 
they can be easily calculated once they are related to the inviscid 
pressure by setting 0 = 9^aU_ tangent-wedge pressure law. 

However, it was found that the variation in pressure across the dis- 
placement body was very large for the compression ramp study, and it 
was reasoned that some advantage might be gained by including the 
higher- ordered terms in the last expression in equation (3.2-6). For 
this reason, it was decided to take 

(If) (3.2-7b) 

\ /s, o 

where pg and Ug are the inviscid surface values which include 
displacement -body effects. 
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Upon substitution from equation (3-2-7b), equation (3- 2-5) 
now be written as 


Pe = PdB - “^PeUe^S . (5.2-8) 

Here and Ug depend implicitly on as will be shown later 

and an iterative solution will be employed to extract the value of 
Pg from equation (3.2-8). 

Note that for the flat -plate study the curvature is zero and 
equation (3.2-8) gives the result that to second order P^ = Ppg . 

In order to relate Tg, Pg, and Ug to the inviscid sur- 

face pressure P^ , it will be assumed that the shock layer is 
isentropic as was done by Flugge-Lotz and Blottner (1962) as well as 
most other investigators of the problem. 

Note that for the present model, the inviscid pressure on the 
displacement body is given by the tangent-wedge pressure law which 
although being approximate is nevertheless generally accepted as being 
valid. In addition to the inviscid pressirre, it is necessary to have 
two other properties with which to define the isentropic flow. The 
second property is 


He = Hoo (3.2-9) 

which in no way limits the flow model since it is exact. 

The remaining property sought is the isentropic constant here 
taken to be the free- stream value given by 
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Pg/Pe = Pco/pI • (3.2-10) 

This assumption is not without precedent since all previous authors 
have adopted it. The alternate approach of evaluating the isentropic 
constant directly behind the leading-edge shock is unacceptable since 
the effective shock angle is not known aproiri. However, it should 
be pointed out that downstream along the body the outer reaches of 
the boundary layer merges with an inviscid flow that has crossed a 
much weaker portion of the leading- edge shock and as such the effective 
isentropic constant in this portion of the flow would be expected to 
be nearer to the free-stream value. Thus, the inviscid surface pres- 
sure is taken to be related to the temperature through the relation 

Pe/P.» = (Te/Too) A . (3.2-11) 

Assuming that Pg/P^^ computed from the tangent-wedge relation, 
equation (3.2-11) may be solved for Tg to obtain 

Izl 

Te = T^(VPoo) ^ ^ (3.2-12) 

which can be used to evaluate pg by either the Sutherland viscosity 
law (equation (1.5-6)) or the linear viscosity law equation (I.5-7). 
Substituting for T^/t^ in equation (3-2-12) from the equation 


of state, we find that 
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1 

Pe = (Pe/Pco)’'. 


(3.2-13) 


where it is noted that = 1 in the present nondimensional scheme. 

To obtain Ug as a function of we note that on the 

body surface 

He = Te + uf/2 (3.2-14) 

from the total enthalpy relation. If we first write the total enthalpy 
constant Hg as 

He = + U^2 (3.2-15) 

and note that = 1 and Too = ^ — n, equation (3.2-14) can be 

(7-1)M^ 

solved for Ug to find that 


Ue - 




1 + 


1 - (Pe/Poo) ^ 


2 




(3.2-16) 


Thus^ we have shown through equations (3.2-12-16) that assuming an 
isentropic flow behind the shock, Tg, Pe. Pe^ and Ue can be related 
to the inviscid surface pressure Pg. 

Similarly dUg/d^ can be related to the pressure gradient by 


first noting that 
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dUg 
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d^ 

d^ 

ds 

d| PgdeUe 

ds 


(3.2-17) 


To relate the pressure gradient dPg/ds to the velocity gradient 
dUe/ds , we observe that the inviscid s-momentum equation when written 
in terms of the present coordinate system is from Van Dyke (1962c) 


ds 


-P 


U — + V 
ds 



(3.2-18) 


Employing the irrotationality condition given in equation (l. 2-l6), 
equation (3.2-18) when evaluated at the wall becomes 


dPe 

ds 


= -Pe 


U. 


dUe 

ds 


+ V, 


/8v' 


(3.2-19) 


Ordering terms in the above equation according to equation 
(3. 2-5), we obtain that to second order 


d^ 

ds 


p dU„ 

^ = -PeUe-7^’ 
^ ds 


(3.2-20) 


which allows equation (3.2-17) to be written in the desired form 


dUg _ -dPg/ds 

d| p2|j ij2 

e e e 

'•'The analysis presented here is for the plane flow case. 

J = 0. 


( 3 . 2 - 21 ) 


Thus 
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Returning to equation (3.2-8), we may now write it in the 
following form once substitution has been made for and from 

equations (3-2-13) and (3-2-16) 



To determine the surface pressure gradient dP^/ds in terms of the 
displacement- body pressure gradient, we differentiate equation (3-2-22) 
with respect to the s-coordinate and solve for dPg/ds to obtain 




( 3 - 2 - 23 ) 

Note that equation (3-2-22) cannot be solved explicitly for the 
inviscid pressure Pg . Thus, when equations (3-2-22-23) are applied 
at a station, an initial guess of Pg/P<„ will be used to evaluate 
equations (3-2-22-23) for the first iteration, thereafter the value 
of Pg/Poo from the previous iteration will be used until the overall 
solution and thus Pg/P^ converges. 

To compute dP^g/ds necessary for evaluating equation (3.2-23), 
the tangent-wedge pressure law is differentiated to give 
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dPoB 7PcoMoc(eS" + 9 ’wall) 


ds 


(1 + b^m2g2)V^ 


1 + 2B^1V^0^+ 2BM„efl + ^ 


where 
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(:$.2-24a) 


= 0 . 


wall 


+ £ 


ds 


(3.2-24b) 


and 


B = 


7+1 


(3-2-24c) 


In summary, the inviscid properties necessary to compute a 


boundary-layer solution at a particular station are Pg, Tg^ f^e^ Pe^ 
dUa 

U_ , — ^ , and I . This section has shown how these are related 

to the inviscid surface pressure Pg and its gradient dPg/ds and 


as such the problem is now reduced to one of the determining Pg and 
dPg/ds. Furthermore, through equations (3*2-1) and (3*2-24) the 
values of Pg and dPg/ds can be seen to depend only upon the unknowns 
6' and 6" . Thus, once 6' and 6" are known all inviscid pro- 

perties can be computed. 

3. 3 Solution at a General Station 

As was noted in Section 3*2, all inviscid properties can be 
calculated once 6' and 5" are known. Therefore, the overall 
method of solution can be s\immarized as follows. (l) Make an initial 


Note that ( )'= — and this notation will be used 
ds 

interchangeably. 



8l 


estimate of 6' and 6" . With these estimates calciilate all the 
inviscid properties and make one iteration on the bound ary- layer solu- 
tion. (2) From that houndary-layer profile , integrate equation 
(1.5-7) and determine the displacement thickness 5.^ Employing 
finite-difference formulas calculate 6* and 5". (5) Repeat steps 

(1) and (2) until the solution has converged within some acceptable 
error. 

a. Initial Estimates of Inviscid Properties . - To apply the solu- 
tion technique at some general station, we will estimate the inviscid 
properties necessary to compute a boundary-layer profile by adopting 
the following estimates on 6' and 5". 

~ ^6s ^OS ^ (5'3-la) 

and 

Sns = 5qs + Sqs ( 3.3-lb) 

where As is the step size in the s-coordinate and the subscripts 
NS and OS refer to the new station and old station respectively. 

b. Interaction Technique .- Using the initial guess at the bo\mdary- 
layer profile described immediately above, equation (1.5-7) is 
integrated to determine the new estimate of the displacement thickness 

6. With the following two -point derivative laws 

"^See Appendix A for details of method employed to calculate S 
once 5 is determined from equation (1.5-7)- 
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— - 52)/Asp, 

(3.3-2a) 

= (6| - 5^)/Asi, 

(3.3-2b) 

= (5"- 5p/As^, 

(3.3-2C) 


or the three-point derivative laws 


Asi(Asi + Z^sp) 


62 ^ 55 

As 2 

_ H As^ (Asp + As2)^/ 

Asp (ASp + As2)^ 


(5. 5- 5a) 


2 Asp(Asp + Asp) 

S2_ _ 53 5^ 1 

' 1 1 \ 

As2 

As^ (Asi + As2)^ ^ 

yAs^ (Asp + As2)^/ 


+ 



1 

(Asp + Asp) 


"/i 


(5. 3- 5b) 


Asi(Asi+ As2 ) 


- + 

] 

Asp 1 

1 

^ \Asp (Asp + Asp)^/ 

Asf 

(Asp + Asp)2 j 
J 


(3. 3- 3c) 

the displacement thickness derivatives are updated. In the above 
difference formulas the error in the two-point law is 0 (As). The 
advantage of the three-point laws is that the error in the 5’ 
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formula is O(As^). Figure A shows the grid for both laws which have 
been written for a variable step size. 


Downstream 



Figure A 


Using the new estimates of S' and S" the inviscid properties can 

be updated and another boundary-layer profile calculated and this 

cycle one woiild hope to repeat until the solution converges, that is, 

stops changing within some criteria. However, it was found that the 

above iteration process was divergent which is in agreement with the 

results of Fliigge-Lotz and Blottner (1962). The technique they 

employed to obtain a convergent scheme is as follows. Assuming a 
+ 

value of Pe/^00 > inviscid properties were computed and a boundary- 
layer profile was calculated from which a new estimate of Pp/P was 
made. Assuming a second value of P^/p^, the above process was 
repeated. Then with the assumed and corresponding calculated values 
of Pg/Poo between two consecutive iterations, a linear extrapolation 


was made to determine the value on an assumed versus calculated plot 

where ' Vf»)ealculatea- ««> this walue 

as the new assumed value, the above cycle was continued tintil the 


Fliigge-Lotz and Blottner (I962) computed the pressiire gradient 
with a two-point difference formula between the pressure at the new 
station and the pressure at the previous station. 
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actual calculated value of Pg/^oo equaled the assumed value within 
some acceptable error. 

The method employed here to obtain a stable interaction process 
is very similar to that en^jloyed by Fliigge-Lotz and Blottner (1962). 
At a new station the initial estimate of S"* is held fixed when 
calculating the inviscid pressure gradient on the displacement body 
(equation (3- 2-24)) until the calculated value of 6" stops changing 
to some arbitrary level which in this study required 


^"/S'ave < 


(3.3-4a) 


where 


= 


6" (present iteration) - 6" (previous iteration) 


(3. 3-4b) 


and 


ave 


^ 1 
2 


S" (present iteration) + 5" (previous iteration) 


( 3 . 3-4c) 

Although the 6" value used to compute dP^g/ds is held fixed in the 
above process, 6' is allowed to adjust according to equation (3.3-2a) 
or (3*3-3a). Thus, the inviscid properties are updated after each 
iteration. 

The value of 6" which was held fixed when computing the pressure 

gradient we will call 6^, the subscript indicating a guessed value. 

The calculated value of 5' corresponding to d" will be called 6" 

G C 
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meaning calculated. If we now make a plot of all versus 5" 

b- U 

points, then it can be seen that the desired solution has been found 
at a particular station when 6" = • This being true, we need only 

0 Cj 

to check to see whether Bq = 6q within some prescribed criteria. If 

^Gc/^GC - (3.3-5a) 

where 

^GC ~ 8" (calculated value) - 5" (guessed value) , (3- 3-5^) 

and 

^GC “ 5 8" (calculated value) + 8" (guessed value) , (3*3-5c) 

is met, the solution was considered to have converged at a particular 
station. If equation (3-3-5) was not satisfied then the present 
^G'^C previous one were used to linearly extrapolate 

what the desired value of 8" sho-uld be. With this extrapolated 
value of 8" as the new 8^ value, a new 8^ value was computed 
and equation (3-3-5) was checked again. The process of extrapolating 
the two previous values of 8^ and 8|^ to obtain a new 8^ value 
was continued until equation (3-3-5) was satisfied.’*' When satisfied 
the program moves on to the next station. Figure 5. 1 shows the varia- 
tion of 8^ and 8|( for the = 9.6 flat-plate case. In all 

'“'This method was found to be stable for accuracies as small as 
10~° in equation, (3-3-5a). 
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examples this variation was found to be approximately linear and the 
most guessed values necessary was 3 oi" 

The station at which the interaction is initiated is special in 
that the above technique cannot be applied there. The methods employed 
to obtain the initial profiles are given in the next section. 

3 . ^ Initial Profiles 

As was noted in Chapter I in order to begin the numerical integra- 
tion of the boundary-layer equations, the solution must be known at 
the initial station. Thus the profiles at the initial station must be 
fed in or computed in some manner. 

The present study will be limited to bodies with a sharp leading 
edge whose leading-edge surface is approximately flat and at zero angle 
of attack relative to the mainstream. Thus, near the leading edge 
where the interaction will be initiated, we will adopt for our analysis 
the boundary-layer equations without curvature effects. Two similarity 
methods will be put forth with which to compute the Initial profiles. 

The first considered is the classical strong interaction solution here 
called the hypersonic limit method, and the second its modified form. 

a. Hypersonic L imit Method . - The similarity method presented here 
for obtaining initial profiles and inviscid properties is basically the 
one given by Lees and Probstein (1955) who obtained a first-order esti- 
mate to the displacement thickness and surface pressure by expanding 
these quantities in asymptotic series and calculating the leading term 
in the expansions. The first-order estimate obtained by Lees and 
Probstein (1953) was subsequently adopted by Flugge-Lotz and Blottner (l962'i 
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in their work with insulated walls. The reason for the present version 
is to allow for a wider variation in wall temperature. 

With this method, we note following the analysis of Hayes and 
Probstein (1959) or Sullivan (1969) that certain simplifications 
become apparent if one considers the case where the inviscid Mach 
niimber Mg tends toward infinity. Writing equation (3.2-14) in terms 
of the inviscid Mach number we obtain that 

rp 

- i > 0 as Mg 00 , (5- 4-la) 

Hg 1 + 2Ll_i 

e 


2T 

e _ 

2 

ur 2 _ 


-i M? 


0 as Mp -> 00 , 


(3. 4- lb; 


and 


2H, 


7-1 ,,2 

Me 


+ 1 1 


as Mg -» <» 


(3. 4-lc) 


If we assume also a linear temperat-ure-viscosity law thus making 
Z = 1 and set the terms equal to zero in equations (1.4-8-10), 

the limit form of the governing equations for the self- similar flow 
are obtained as 
CONTINUITY EQUATION 


— + F = 0 , 
dTi 


(3.4-2) 
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S -MOMENTUM EQUATION 


F" - VF' - p(G - p2) = 0 , 


(5-^-3) 


and the 


ENERGY EQUATION 


G" - oVG' + 2(o - 1)(FF')’ = 0, 


(3.4-4) 


where 3 which is given by 


3=1 + 


2Tg 


U„ d? \2T, 


2 '• 


(3.4-5a) 


in the hypersonic limit goes to 


3 



/ — dx ^ . 

- ° (Pe/P~)^ 


( 3 . 4-5b) 


This last resiilt is obtained by first employing the definition of 5 
along with the hypersonic limit value of Ug ~ constant implied by 
equation (3-4-lc). Thus, it can be seen that in order to obtain the 
initial profiles from a numerical solution of equations (3* 4-2-4), 
one must first determine the proper value of 3 in equation (3- 4-3)* 
Following the lead of Shen (1952) and Hayes and Probstein (1959) ^ 
we adopt the following pressure model 



s-1/2 




where C is a constant which must be determined. 


(3.4-6) 
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Substitution of equation (3*^-6) into equation (3*^-5ti) leads to 
the familiar result ('see for example Hayes and Probstein(l959) ) • 




(5.'^-7) 


Thus, with the value of (3 given by equation (3*^-7)^ it can be 
seen that solutions to equations (3* ^-2-4) depend only upon the 
boundary condition 1 1 = which when supplied allows one to 

compute the initial profiles. 

In addition to the initial profiles, one must also know all the 
inviscid properties at the initial station, all of which can be com- 
puted from equation (3- ^-6) once C has been determined. 

To compute C , we first observe that when curvature effects are 
neglected, equation (1.5-6) may be written as 


S = 


.111 f“j'x . iL' a,,, 

o \Te < 


u 


e '^e 


(3.4-8) 


which becomes in the hypersonic limit 


= ^ ~ mJ \! 


2M-cx,Co 


P /P dx 
o e' oo 


1/2 


2 P /p 

e' 00 


where 


Q = (7 - 1) m' 




3/4 


00 y — *5, S J 


Q = / (G - p2) dq 
^ O 


(3.4-9) 


(3.4-10) 
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The derivative of the displacement thickness given in equation (3-^-9)> 

is now substituted into the tangent-wedge pressure law and the resulting 

pressure set equal to that obtained by equation (3*^-6) to give in 

the limit of M — > OO 
00 


C s 


- 1/2 



(3-h-ii) 


where near the leading edge we have made use of the fact that 


= , ( 3 . 4 - 12 ) 

and B is given by equation (3’2-2kc). After substitution for 6 
in equation (3- ^-11) from equation (^.h-9), we find that 


= 7 (7 - 
4 



( 27 B)^/^Q 


1/2 


(3.4-13) 


where 



p*U*s* 
00 00 


* 




00 


7 


(3.4-14) 


which now allows us to write equation (3-4-6) in the standard form 

Pg/P,, = a X , 


(3.4-15) 
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where here 


and 


a = I (7 - 1) y^ 27 B Q,, 


(5. ^-16) 



{3-^-lD 


The only inviscid property at the initial station not known but 
required for later solutions is 5 which can now be computed by 
substituting equation (3* ^-6) into the definition of integrating 
to obtain 


I 


- 2p^C 



(5A-I8) 


In equation (3- ^-15)^ k is recognized as the hypersonic interaction 
parameter which indicates large displacement induced pressures in 
regions of large X . Thus, regions of large X have come to be 
called strong interaction regions whereas regions of small X have 
come to be called weak interaction regions. Note that once the wall 
temperature is specified, equations (3* ^-2-U) can be solved numerically 
and Q (equation (3.4-10)) can be evaluated thus deterraing ct for 
that particular wall temperature. Hayes and Probstein (1959) present 
tables of a versus wall temperature for a representative range of 
flow variables as obtained from an analysis similar to the above but 
originally presented by Lees and Probstein (1953). 
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b. Modified H^ersonic Limit Method . - In the discussion of results, 
it will be shown that a blimp appears in the computed surface pressure 
gradient when the above hypersonic limit method is employed to obtain 
the initial profiles and inviscid properties. Because the bump 
appeared immediately downstream of the initial station, it was reasoned 
that one source of it was due to the fact that at the initial station 
the inviscid properties were the hypersonic limit ones whereas at the 
succeeding stations this limit is not applied. That is to say, when 
computing the inviscid properties at the second station, there are 
terms included which were considered of higher order at the initial 
station where the hypersonic limit was applied. Therefore, in order 
to remove one possible source of the bump, the modified hypersonic 
limit method was formulated. This formulation consists of keeping 
higher- ordered terms when computing the inviscid properties at the 
initial station, thus making that station more consistent with the 
downstream stations. As before, we will take the inviscid surface 
pressure to be correctly given by equations (3-^-15) and (5.4-16-17) 
where Q in the modified method is given by 


Q = 




dT) , 


(3.4-19) 


which in the hypersonic limit can be seen to reduce to equation 
(3.4-10). In addition, we will assume also that the governing equa- 
tions for initial profiles are correctly given by the limit equations 
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(^.k-2-k) and ( 3 .^- 7 ). Thus, it is seen that the initial profiles 
for both hypersonic limit methods are identical and the only 
difference is in the method employed to obtain the inviscid properties. 

To compute a from equations (3.^-l6) and ( 3 . ^-19)^ ^ iterative 
approach will be taken. Assuming the initial profiles have been 
mimerically obtained, a is estimated and the necessary values of 
Tg and Ug to numerically integrate equation ( 5 * ^- 19 ) are determined. 
With this value of Q, a new estimate of a is made and the cycle 
repeated until a is calculated within some acceptable error. 

Once a is determined, all inviscid properties with the excep- 
tion of 5 are known. To compute | in the modified method, we 
will set the relation of p given in equation (5.4-5a) equal to its 
hypersonic limit value and solve for | . Thus, 


5 = - 



2 

^e^e‘ 


e e 


I 1 + 


dPe ^ 
Ue2y ds 


(3.4-20) 


which can now be computed. The value of 6 necessary to compute S' 
at the next station is obtained from equation (5.4-8) whereas S' is 
determined by equating equation (5*4-15) to the tangent-wedge pressure 
law and solving for S' . 

To compute the additional values of S and S' necessary to 
apply the three-point difference laws at the next station the 
following procedure was adopted. For the modified method S' was computed 
as stated except s is given a value two As steps upstream of the 
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desired station. The additional value of 6 for the modified method is 
obtained from a two-term upstream Taylor series expansion about the 
initial station. For the hypersonic limit method S and its deriva- 
tives were computed from equation (3.4-9) by taking a value of s two 
As steps upstream from the desired station. 

3> 5 Discussion of Results 

a. Flat-Plate Studies . - The first case investigated was that of 
flow over a flat plate with displacement thickness interaction. This 
case was originally investigated employing finite-difference methods by 
Fliigge-Lotz and Blottner (1962). Thus, this case will serve as a 
test case with which to compare the present study. 

The other justification for studying this case lies in a con- 
troversy over whether or not branching solutions result when problems 
of this nature — ^particularly the flat-plate case — are investi- 
gated numerically. Garvine (1968) observed a saddle-point (or 
branching) behavior with his viscous shock-layer solution for the 
flat plate. Similar behavior was noted by Baimn (1968) with his super- 
sonic boundary-layer solutions. However, Flligge-Lotz and Blottner 
(1962), Siillivan (1969), Stollery (1970), and Cheng, Chen, Mobley, 
and Huber (1970) did not observe any tendency toward branching solu- 
tions and thus a reworking of the hypersonic flat-plate case seems 
appropriate as a test of the generality of Fliigge-Lotz and Blottner 's 
technique. 
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1 . Initialization Procedures and Effects . - In Section 3- two 
methods were given with which to compute the initial profiles neces- 
sary to begin the interaction. The method adopted in this study was 
the second of these - the modified hypersonic limit method. The 
reason for selecting this method can be seen from Figure 3* 2 where 
the second derivative of displacement thickness - the most sensitive 
solution function, is shown for the case of a plate in a Mach 12 free- 
stream. Figure 3-2 shows an oscillatory type behavior developing when 
the hj^ersonic limit method is employed to obtain the initial profiles 
but shows a smooth transition when the above modified scheme is 
employed. Since the pressure is a function of 6' it becomes apparent 
that 6" reflects the behavior of the pressure gradient in the 
s-direction. Thus, to assure a smooth pressure gradient in the 
s-direction for the first several stations the modified method was 
adopted. Figure 3* 3 shows the wall pressures corresponding to the 
above case for both hypersonic limit methods. The differences in 
wall pressure between the two methods can be seen to fall outside of 
the plotting resolution implying the oscillations of Figure 3. 2 are 
apparently of only secondary importance to the overall solution. 

Figure 3.4 gives an indication of the effect of step size on the 
solution near the initial station in terms of how it influences 6" . 
From equation (3-4-9), we obtain upon differentiation that for the 
strong interaction region (large X) 


6 " 


■5/4 




(5.5-1) 
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Thus, on Figure 3.4 equation ( 3 . 5 -I) should appear as a straight line 
and with this in mind step sizes of 0.0025 and O.OO 5 were chosen 
since they yielded the approximate desired behavior. 

Figure 3- 5 shows the step size effect on the computed surface 
pressure to be very small and to diminish a's the solution proceeds 
downstream again pointing to the basic consistency of the present method. 

2. Downstream Instabilities .- As the solution marched downstream 
there appeared an oscillatory behavior in the computed surface pressure 
gradient for the flat- plate case (and the compression ramp case dis- 
cussed later) resulting in an unwanted termination of the computation. 

A similar behavior was observed by Cheng, Chen, Mobley, and Huber 
( 1970 ) with their fully viscous shock-layer solutions in which they 
found the instability to enter through the continuity equation. They 
reasoned that its cause was due to the diminishing wall slip influence 
allowed in their solution which resulted in a loss in accuracy in the 
technique used to integrate the continuity equation. The instability 
in this study also seemed to originate with the continuity equation 
although no concrete explanation can be offered as to its origin. 

Initially, it was believed that the oscillatory behavior was 
associated with the method employed to obtain the initial profiles 
and/or related to the original step size employed. Therefore, the 
first modifications made were those reported in the previous section. 

When those modifications failed to eliminate the instability it was 
then felt that the oscillatory behavior might be related to errors 
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present in the finite-difference derivative laws used to compute 6' 
and 6" . Thus, a study was made to determine if the more accurate 
three-point derivative laws would be more appropriate. Figure 3- 6 
shows the computed second derivative of displacement thickness 
employing both laws and it can be seen that the three-point difference 
laws allows the solution to proceed further downstream to an s of 
0.11 before running into the stability problem again. From Figure 
3. 6, it became apparent that although the derivative difference scheme 
exerted some influence on controlling the instability and apparently 
eliminated it for a while, another method would have to be found to 
abolish it from the problem. 

Taking a hint from the identification by Chen et al. (1970) of 
the continuity equation as the source of the instability, the method 
adopted here for filtering out this oscillatory behavior consisted of 
simply averaging the V-velocity profile - as it was being calculated - 
with its counterpart for the same value of one station back along 
the body. Such a technique is formally valid to the same order as 
the difference equations, that is, to order As and therefore intro- 
duces no new errors to the solution. Figure 3.7 shows the computed 
surface pressure for the M<x, = 9-6 flat-plate case for both the 
solution with and without the averaging technique. The pressirre 
gradient behavior can be seen in Figirre 3. 8 where the second deriva- 
tive of displacement thickness is shown to vary smoothly over the 
entire region considered for the filtered solution. 
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Since the Mco = 12 flat-plate case showed no oscillatory behavior 
this case was run with and without the averaging technique to ascertain 
what effect the averaging had on the computed solution. As might 
have been expected, the differences between the two solutions was 
virtually undetectable. 

5. Comparison of Experimental and Theoreticl Flat- Plate ResilLts 

Example 1 

The first case studied was the Moo = 9*6 flat- plate case pre- 
viously investigated by Flugge-Lotz and Blottner (1962). The flow 
conditions for this case are given as 


M =9.6 
00 

(3. 5-2a) 

T* = A2. 
“00 

( 5 - 5 - 2 b) 

II 

H 

(4.5-2c) 

a = 1.0 

(3. 5-2d) 

, /s* = 0.846 X 10 ^ 

'00 

( 3 . 5-2e) 

T /T = 1.0 
w 0 

( 3 . 5-2f) 

L*= 1 inch. 

(3.5-2g) 


Figure 3.9 shows the pressure induced on the plate surface when a 
linear viscosity-temperature law is employed to compute the numerical 
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solution. Also shown are favorable conqjarisons with the results 
reported by Flugge-Lotz and Blottner (1962) and the experimental data 
of Bertram and Blackstock (196I). The present solution was initiated 
at s = 0.01 which corresponds to a X of approximately 25 and is 
therefore in the strong interaction regime. Figure 5-10 shows the 
above displacement thickness to be slightly lower than that predicted 
by classical boundary-layer theory. As s becomes large, the skin 
friction behavior can be seen from Figure 3.11 to approach that 
predicted by classical boundary-layer theory. 

Example 2 

The flow conditions for the second example are 


M = 12 
00 

( 5 . 5-5a) 

T* = 596o°r 

(5.5-5b) 

7 = 1 . 4 

(3.5-5C) 

a =0.70 

(5-5-5d) 

Re- /S* = 0.79 X 10^ 
■^00 

(3.5-3e) 

TjTo = 0.15 

(3.5-5f) 


and 


= 1 inch . 


(^•5-5g) 
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This case was previously investigated by Sullivan ( 1969 ) who obtained 
an approximate solution by employing Lees ( 1956 ) cold- wall similarity 
analysis. The interaction for the present example was initiated at 
s =0.01 which corresponds to a X of approximately 58 * 

Figiire 3.12 shows the pressure induced on the plate surface when 
a linear temperature-viscosity law is employed for the numerical 
solution. Also shown are the approximate solution of Sullivan ( 1969 )^ 
the second- order weak interaction theory reported in Hayes and 
Probstein (1959) ^ and some experimental data by Hall and Golian (i 960 ) 
along with two points from the present study employing Sutherland’s 
viscosity law-equation (l. 3-6). As can be seen from Figure 3*12, 
for X < 6. 5 the present solution employing a linear ten^erature- 
viscosity law virtually reproduces the second- order weak interaction 
result. The differences between the values of wall pressure obtained 
by employing the linear law and Sutherland's viscosity law indicate 
only a small influence due to the viscosity law employed. 

Shown in Figure 3*13 is the heat transfer to the wall along 
with the experimental data of Hall and Golian (i 960 ). The agreement 
appears acceptable. 

The displacement thickness and skin friction for the above case 
are shown in Figures 3- 1^-15 and for large s are observed to 
approach the behavior predicted by classical boundary- layer theory 
thus indicating a diminishing displacement- body effect. 
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b. Compression Ramp Studies 

1. Geometry and Flow Conditions.- The locus of points on the 
body surface are given by y = x5/l50 up to x = 5-1 inches where 
the cubic surface merges with a 28 ° wedge. 

The flow conditions for this case are taken from Stollery (1970)^ 
who in addition to experimentally investigating the flow also calculated 
some approximate solutions with various approximate inviscid pressure 
laws. However, in his analysis, curvature effects were neglected and 
one of the purposes of this chapter is to determine whether curvatiire 
effects are important for problems of this type. The flow conditions 
studied here were 


Moo = 6 , 8 , and 12.25 

( 3 . 5 -^a) 

T* = 2340°R 

( 3 . 

1 — 1 
11 

(3.4-4c) 


a = 0.72 ( 3 . 5 -^d) 

Re„ /s* = 0.858 X 1q5, 0.458 X 105, and O .258 x 10^ (3- 5-4e) 

^co 

i;^/To= 0.223 (3.5-4f) 


and 


L* = 1 inch . 


(3-5-4g) 


Sutherland's viscosity law was employed for the ramp study. 
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2. Displacement Effects. - The first study was made at Moo = 12.25 
and was aimed at determining the effects of the displacement body on 
the wall pressure for the case where curvature effects were neglected^ 
hereafter referred to as the displacement problem. Figure 3.16 shows 
the wall pressure that results from flow over the displacement -body 
when curvature effects are neglected. Also shown in Figure 3*l6 is 
the wall pressure when displacement -body effects are ignored, that is, 

0 = Q^all ‘ seen from Figure 3*l6>the wall pressure which 

includes displacement- body effects predicts the trend in the experi- 
mental data very well up to an x of approximately 2. 5* From x ~ 3 
on, the displacement solution approaches the wall pressure where dis- 
placement effects are ignored thus indicating that the local body 
®wall dominating the displacement- body effect. This 
behavior is verified by the displacement thickness distribution of 
Figure 3> 17 ~ showing a marked thinning of the boundary layer aft of 
the X = 3 point. Figure 3.17 also shows a comparison between the 
displacement thickness obtained in the present study by considering 
displacement- body effects alone and the approximate theories of 
Cheng et al. (196I) and Siillivan (1969) as applied by Stollery (197O). 
Thus, it is encouraging to find the present solution closely following 
that given by Stollery using Sullivan's method and clearly delineating 
a marked thinning downstream on the body. Figure 3.18 gives the wall 
heat transfer ratio 9/q,F.p, where q is the heat transfer predicted 
by classical bound ary- layer theory which was found for this case to be 
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qp.p. =0.2503 , (3«5-5) 

along with the experimental data ofStollery (1970). The solution which 
includes displacement-body effects appears to predict the trent 
reasonably well up to x = 2 but thereafter drops below the experi- 
mental results. 

As was mentioned previously in Section 3-^s,(2), the instability 
discussed there appeared for the above case when the averaging tech- 
nique on the V-velocity profile was not employed. However, in this 
case the instability appeared at a much larger value of s(x ~ 1.86) 
and the point at which it appeared was found to depend on the step 
size As . For example, when the step size was doubled aft of x =1.1 
over the previous step size, the instability did not appear until 
X 2. 32. Thus it appears that a second possible method of elimina- 
ting the instability woiild result from a careful selection of the 
step size increment (As). 

3* Curvature Effects . - In the study of the effects of surface 
curvature on flow up the Mach 12.25 compression ramp, an attempt was 
made to analyze curvature as it affects the inviscid properties by 
allowing for a pressure variation across the displacement body, 
independent of curvature as it affects the viscous region by also 
allowing for a pressure variation across the boundary layer. 

The first such investigation sought to determine the effect on 
the interaction problem of allowing for a pressure variation across 



104 


the displacement body with no variation in the bound ary- layer pressure. 
Figure 5-19 shows the wall pressure that results from that case where 
only the inviscid curvature corrections associated with the displace- 
ment body were considered, that is, the pressure across the boundary 
layer was taken to be constant while the pressure across the displace- 
ment body was taken to vary according to equation (3-2-22). Also 
shown in Figure 3*19 is the wall pressure for the displacement problem 
alone. From Figure 3*19^ it can be seen that the curvature correction 
is surprisingly large - Figure 3*20 showing the increase to be on the 
order of 20 to 50 percent in the region of interest. Figure 3*20 shows 
that the pressure correction closely follows the local variation of the 
ratio of displacement thickness to the radius of curvature (or kS). 

It is quite clear then that even with the relatively small values of 
K 6 shown in Figure 3' 20, the inviscid normal pressure gradient 
encountered implies sizable normal pressure gradients are to be 
expected in the viscous region therefore indicating our solutions are 
necessarily incomplete at this point. 

When an attempt was made to analyze boundary-layer curvature 
effects without first adjusting the inviscid flow properties for 
displacement or curvature effects (i.e., we took 6 = 9^ ^ ) , a 
branching type behavior was observed in the flow properties as the 
solution marched away from the leading edge. Figure 3.21 shows the 
above resulting wall pressiire distribution dropping off towards 
zero (corresponding to an expanding flow) when the computation 
terminated. The cause for this behavior seems to be the low 
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estimates of xnvlscid surface pressure near the leading edge obtained 
by taking 0 = 6-^a.ll coupled with the increasing variation in pres- 
sure across the boundary layer - approximately 13 percent at the 
point of termination. Noting that the surface curvature is negative, 
it can be seen from equations (l. 2-ii-) and (l. 2-6b) that the maximum 
pressure in the boundary layer must occiir at the wall which from 
equations (l. 3-8g), (l. 3-8e), (l. 3-8i), and (1.4-11) is obtained as 

Pwall = Pe + Ugl^all ‘ (3.5-6) 

From equation (3.5-6), it is observed that if the estimate of is 

low^ then there will be a correspondingly low value of wall pressure 
in the boundary layer. When this is coupled with a large variation in 
pressure across the boundary layer, the resulting effect is to cause 
low, unrealistic densities to occur in the outer reaches of the viscous 
region. This in turn results in a subsequent increase in displacement 
thickness that would have been offset had interaction effects been 
included simultaneously. However, since displacement effects are not 
included there exists no mechanism by which the boundary layer can 
adjust to the real problem. We note that a method was found to filter 
out this expansion branch and consists of averaging the pressure 
integral in equation (l.4-l6) in exactly the same manner as the 
V-velocity profile was averaged to eliminate the oscillatory 

'*'In the region of concern. Figure 3«l6 shows this estimate 
(0 = ©wall) 50 to 90 percent lower than that predicted by the 

displacement problem solution. 


instability. Figure 3*21 shows the wall pressure obtained when the 
averaging technique is employed along with the inviscid surface 
pressure obtained as stated by taking 0 = ®vrall ‘ Exactly why this 
correction works is unknown but it is worth repeating that the inclu- 
sion of displacement effects would itself relieve the problem and is 
therefore not of immediate concern to the overall effort. 

When displacement- body effects were included for the above case 
in addition to the boundary-layer curvature effects already accounted 
for^ the branching type behavior which occurred near the leading edge 
was overcome temporarily but reappeared again at approximately 
X » 1.57* The reason for its reappearance is due once again to low 
estimates of the inviscid surface pressure. However, for this case 
the cause of the low surface estimates of Pg is due to the large 
variation in pressure across the displacement body which exists but 
has not been accounted for here. Figure 5.20 indicates the inviscid 
surface value of pressure at x = 1. 5 should be approximately 50 per- 
cent higher than the displacement-body value which was used for this 
case. With the inclusion of this final surface curvature effect, we 
come to what will be referred to as the complete second~order solution. 

2. The Complete Second -Order Solution . - When all the second- 
order effects due to displacement body and sirrface curvature were 
accounted for, no branching type behavior was observed in the flow 
properties. Figure 5-22 gives the wall pressure for the complete 
second-order solution along with the wall pressure predicted by 
classical or first-order boundary-layer theory, that is, neglecting 
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displacement and curvature effects. Near the leading edge where the 
surface curvature is small. Figure 5*22 shows the displacement effect 
to be large and to diminish as the solution proceeds downstream. In 
addition. Figure 22 implies that although the effects of displace- 
ment speed and longitudinal curvature cause* large pressure gradients 
to exist across both the displacement body and the boundary layer all 
along the body surface, their effect on the wall pressure is only 
noticeable in regions where the displacement speed is clearly the 
dominant effect. Figure 3*23 shows the trend in the wall heat trans- 
fer data to be predicted reasonably well by the complete solution up to 
X 2. The skin friction for the complete solution is given in 
Figure 3*24 along with the values predicted by classical boundary- 
layer theory and the present interaction technique when curvature 
effects are neglected. Even though we realized that the displacement 
effect was dominant near the leading edge, it was nevertheless quite 
surprising to learn from Figure 3-23 that the curvature effect does 
not noticeably Influence the skin friction distribution for values 
of X less than approximately 2. 5. 

One final test of the technique developed here for studying 
displacement and curvature effects consists of varying the flow 
Mach n-umber and free- stream unit Reynolds number to determine their 
influence on the con^juted wall pressure and heat transfer. Figures 
3.25 and 3-26 show the effect of free-stream Mach number on the com- 
puted wall pressure and heat transfer respectively and it is observed 
that a decrease in Mach number results in a decrease in the computed 
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wall pressure but an increase in the wall heat transfer rate. The effect 
of varying the free- stream Reynolds number while holding the flow Mach 
number fixed can be seen in Figures 3*27 and 3-28 where the corresponding 
distributions of wall pressure and heat transfer are shown. From Fig- 
ures 3-27 and 3*28, it is observed that small decreases in the free- 
stream unit Reynolds number give rise to an increase in the computed 
wall pressure and a reduction in the heat transfer rate. Of spe- 
cial interest is the solution shown in Figures 3- 27 and 3- 28 which termi- 
ated due to separation of the viscous flow from the ramp surface. Figure 
3- 29, which gives the skin friction distribution in the region where it 
goes to zero, shows no unusual behavior as the point of zero shear stress 
is approached - in fact the computation proceeded downstream of the 
zero shear stress point (C^ = 0) before the ntunerical method employed 
to solve the bo\mdary-layer equations became completely unstable due 
to the reversed flow region near the wall occurring downstream of the. 
separation point. In noninteracting boundary-layer theory, it is 
generally accepted that the skin friction derivative dCp/ds has 
a singularity at the point of separation and several studies have 
been made on this phenomena (see for example Brown and Stewartson 
(1969))* There was some doubt as to whether the compressible 
nonadiabatic wall case was singular at separation. However, recent 
studies by Werle and Senechal (1971) have implied that this case is 
also singular. It has been suggested that if the boundary layer were 

\fe note here that since a reversed flow region occurs aft of 
the zero skin friction point, there arises some question as to whether 
the solution is valid beyond that point. 



109 


allowed to interact with the mainstream as done in this study, there 
would not occur a singilLarity in the skin friction derivative at the 
separation point as obtained in noninteracting theory. This behavior 
in the presence of interaction is confirmed in Figure 3- 29 where it 
is seen that no sing\ilarity is observed at the point of separation. 

Thus this verification stands out as an irr 5 )ortant result obtained from 
this study. 

3- 7 Conclusions 

From this study of displacement and curvature effects, the 
following conclusions were reached. 

1. Wo branching behavior was observed for flat- plate interaction 
solutions. 

2. Instabilities may occur in the numerical technique downstream 
of the leading edge unless the continuity equation solution is damped 
by a method similar to the averaging technique employed here. 

3. Accounting for bound ary- layer curvature effects requires a 
good estimate of the inviscid surface pressure in order to avoid low 
unrealistic pressures in hs^ersonic boundary layers. Large inviscid 
and viscid pressure variations were found to exist across the displace- 
ment body and boundary layer for hypersonic flow over curved surfaces. 

4. When the boiindary layer is allowed to interact with the main- 
stream, there appears no singularity in the skin friction derivative 
at the point where the boundary layer separates from the body surface. 
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5> The technique developed here to handle the complete second- 
order theory predicts the trend in the wan properties for regions 
of strong interaction with the mainstream and merges with the behavior 
predicted by classical boundary-layer theory for regions of weak 


interaction. 



IV. GENERAL CONCLUSIONS 


From this study of displacement and curvature effects, the 
fo3J_owing conclusions were reached: 

1. Finite-difference methods can be applied to the composite 
second-order boundary-layer equations to obtain numerical solutions. 

No new difficulties were encountered from such an application that are 
not also present in similar numerical solutions of the first-order or 
classical equations. 

2. For cases where the influence of displacement speed and curva- 
ture are both large, these effects must be accounted for simultaneously 
in order to avoid sizable errors in the wall properties as well as 
numerical difficulties. 

3 . If one is only interested in wall properties, good estimates 
of these quantities can be obtained from numerical solutions of the 
first-order equations taking the inviscid surface pressure to be 
approximately equal to the observed wall pressure. 

4. Comparisons of the composite second- order solutions with 
experiment are respectable. 


Ill 
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VI. APPENDIX A 


NUMERICAL ANALYSIS SECTION 

The method presented here for obtaining numerical solutions to the 
boundary-layer equations is basically the one developed by Fliigge-Lotz 
and Blottner (1962) for solving parabolic differential equations vd.th 
some later modifications by Davis (1970)- The presentation here 
follows closely that of Davis (1970)* 

The s-momentum and energy equations can be written in the general 
form of a parabolic differential equation which is 

^ + ap — + + (x-2 + — = 0 ^ (l-A) 

where i and t) are the similarity coordinates. 

We now replace the derivatives in equation (l-A) with finite- 
difference quotients by first expanding for the unknown quantity, w, 
in the rj direction about the point m, n (see Figure l-A) by means of 
a Taylor series expansion, and then determining the derivatives to 
give a truncation error on the order of the step size squared. 
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Adopting a variable grid step size in the t] direction to allow 
for more grid points in the region near the body where the changes 
are the most rapid, the derivatives at the point m, n for a central 
difference scheme are obtained as 

(VSn)^n^n = [^%-lM%(^% + ^%-l)]wm,n+l 

+ [(At1j^ - ^%-l)/^%^Vl]Vn 

I ^%/^%-l(^% '^%-l)| ^m, n-1 

- i(S\/8T]^) Ati Ai] (2-A) 

6 m,n n n-1 , ' 

and 

(^V^^^)m,n = [2/ATin(Arin + Ari^.i) j w^^ 

- (2/ATinA%.i)wjn,n + [ 2/AT]n_i(ATi 

^ ^%-lfm,n-l - %Ati^-1 

. i(aVaTl^^^^(An^ - ati^_^)2 . 

(3-A) 

The derivative in the i direction is given by a two-point difference 
scheme which has a truncation error of 0(A|) if equation (l-A) is 
evaluated at station m (assuming stations m - 1, m - 2, . . . to be 
known) and 0(A| ) if equation (l-A) is evaluated at the point 


m - 1/2 (Crank-Nicholson scheme). 
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Substitution of these difference quotients into equation (l-A) 
yields a difference equation of the form 

Mn-1+ BnWn + Cn^n-i = (^"A) 

which can be readily solved (see Richtmyer ( 1957 )) once the boundary 
conditions are specified. 

SOLUTION PROCEDURE FOR INCOMPRESSIBLE STUDY 

The overall method of solution as applied to the complete second- 
order incompressible equations is as follows. Beginning at the stag- 
nation point the Sw / term is equal to zero thus reducing equation 
(l-A) to an ordinary differential equation. After making initial guesses 
of all flow profiles and carefully evaluating all invlscid properties 
to avoid indeterminancies, integrate the s-momentum equation using the 
numerical method. Next integrate the continuity and the pressure 
integral equations employing a trapozoidal rule. With these new 
estimates of the pressure integral and velocity profiles update the 
inviscid properties and remaining flow variables and resolve the 
s-momentum equation once again and continue this cycle until the solu- 
tion stops changing within some acceptable error. After convergence 
has been reached step down the body one station and repeat the 
above process using the previous station solution as the initial guess 
at the new station flow profiles. At downstream stations, the Sw / 
is evaluated using the previous station flow profile. 
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SOLUTION PROCEDURE FOR COMPRESSIBLE STUDY 

The overall method of solution for the interaction study differs 
only slightly from the incompressible case. For the complete second- 
order solution we begin at the selected initial station where the 
Sw/^l term is set equal to zero. After m&king initial guesses at 
the inviscid properties and flow profiles, we integrate the energy 
equation using the n\imerical method. With this new estimate of the 
total enthalpy profile; update all viscous properties related to 
temperature. Next, upon solving the s-momentum equation, integrate 
the continuity equation, the pressure integral, and the displacement 
thickness integrals employing trapozoidal rules. After obtaining 
new estimates of displacement thickness and its derivatives, update 
the inviscid properties and return to integrate the energy equation 
thus completing the cycle. When convergence has been reached, the 
program steps down the body one station and continues this process. 

The method employed here to extract the displacement thickness 
from equation (I.5-7) is as follows. We first integrate both sides 
of equation (I.5-7). Now, we must determine at what point the left- 
hand integral is equal to the con^lete right-hand integral. This is 
done by a single Fortran "IF" statement. If the desired value of the 
left-hand integral falls exactly on a grid point, then 5 is simply 
the value of q at that grid point. Should the desired left-hand 
integral value fall between two grid point, a linear extrapolation was 
employed to determine the correct value of displacement thickness. 
Knowing, 6,& was then computed from equation (I.5-8). 
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Step size for the incompressible study 

For both the incompressible and compressible study, a constant 
step size was taken in the t] direction and was Ari = 0.1. For 
the incompressible investigation a step size of As - it/720 was 
chosen along the body surface. 

Step size for the compressible study 

a variable step size was chosen for the compression ramp inter- 
action investigation. For x < 0. 5^ As = 0.0025, for 0. 5 < x < 0. 8 , 
As = 0.005, and for x > 0. 8 was 0.01. The run time for the com- 

plete second-order solution at M = 12.25, RSg = 0 . 858 x 10^ with 
these step sizes was 59 minutes on a GDC 6000 series computer. For 
the M = 6,8 cases, the run-time was approximately 25 minutes. 



VII. APPENDIX B 


THE BLASIUS SHEAR STRESS FOR A GIVEN WALL 
PRESSURE DISTRIBUTION 


The form of the wall pressure distribution is that given by 
Dimopoulos and Hanratty (1968) 

= t 

2 00 00 

where 0 is measured in degrees. 

From Rosenhead (1963), the inviscid velocity is assumed of the 

form 

= U* U^ s^l + bg^ + bi^s^^ + bgs^ + . . . j , ( 2-B) 

The odd terras in equation (2-B) for the assumed symmetrical flow are 
zero as noted by Rosenhead (1963)- The constants U^, b2, b4 . . . 
are unknowns and must be determined in terras of the constants in 
equation (l-B). 

TheBlasius shear stress is given by 
'^^wall 

P*U^2u,2 

00 00 o 

(3-B) 


K-/ ■- 


p*U*U' d*| 

CO 00 o 


= s|E^ + "^Eo + s^Ei, 


+ ^E^ + s®E 


'8 
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where 

E = f"(0), 
o o 


(4a-B) 


Eg = bgfgCO) , (4b-B) 

Ei^ = b^fi;(0) + bgf22(0) . 


Eg = bgf^(O) + b4bgfi;2(0) + bgf^ggCo) , and (4d-B) 

Eg = bgf^(O) + bgbgfg2(0) + b^fi;j^(0) + bj^b^fllggCo) + bgfUi^i^l^Co) . 

(4e-B) 

The fimctions "f" are tabulated in Rosenhead (1963)- We wish 
now to determine the constants bg^ bj^^ , . . . in terms of the wall 
pressure constants given in equation (l-B). 

Noting that for a cylinder 

s* = 0 — d* , (5-B) 

560 


then defining 


360 

gives 



"'"This term was not given by Rosenhead ( 1963 ). 


(6-b) 


(7-B) 
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The inviscid flow is assumed to be governed by Bernoiilli's 
+ 

equation which can be written as 


dP* dU* 

-f = -P*U| — ^ 

d|2$ “ diZ5 


(8-B) 


Substituting the expressions for U* and P* (p* p* ) into 

e e e wall 

equation (8-B), we obtain the following relations for the constants 
4*2^ ^4^ • • • j 


b - 1 ^4 

2 2 2 ^ 
apS 


(9a-B) 




(9b-B) 


^6= ± 


1 ag 

2 a-2 


1 

4 a? 
2 


1/e' 

16 a 5 r 


(9c-B) 


and 


b8 = 


1 ^ 

11 

4 ) 

4 2 " 

_5_ ®'4 5 ^4^6 

2 3-2 

' & 

< a| 

a^J 

128 ^4 16 ^5 




2.1 

2 ^2 J 


/e®. (9d-B) 


The constant Uq given in equation (j-B) is given by 


Here we take uf + ^ uf 

6 0 G 
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Uo 



e 


(10-B) 


We need now to express the shear stress given by equation (3-B) in 
terms of the skin friction used in this study which is 


This works out to be 


_ T* 

C_ = _ wall 


i Pi U*' 


Y 


Re 


(11-B) 


Cf = V 2 


s (Eq + s^E 2 + s^E^ +^Eg + s^E 3 


(12-B) 



TABLE I. - NUMERICAL VALUES OF dUg/ds AT THE STAGNATION POINT 





% 

= dUg/ds 




4o 

174 

175a 

175^ 

177 

200 

200 

Potential theory value 

2.000 

2.000 

1 

2.000 

2.000 

1 

2.000 

2.000 

2.000 

Initial guess equatior 

(2.5-55) 

1.795 

1.654 

1.808 

1.824 

H 

CD 

H 

00 

2.046 

1.57^ 

Final value iterative 
solution 

1.799 

1.645 

( 

1.815 

— 

1.825 

2.042 

1.587 


^Here d*/h* = O.O5 
%ere d*/h* = 0.10 

'''Obtained using c\irvef it A given in Figure 2. Je 
+ + 

Obtained using curvefit B given in Figure 2. Je 
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TABLE II.- COEFFICIENTS FOR EQUATION (2-5-5) 


Red 

4o 

174 

175a 

175^. 

177 

200"^ 

200^"^ 

Co 

-0.558049 


-0. 455695 

-0. 555618- 

-0.408055 

-0. 689685 

-0. 671619 

C2 

-0.766625 


-0. 508756 

-0. 818145 

-0.480901 

-1. 015466 

-0.785560 

C4 

-0. 226754 

-0.441596 

-0. 051854 

-0. 555295 

-0. 051616 

-0. 588269 

-0.048946 

C6 

-0.149195 

-0. 155810 

-0.055420 

-0.101078 

-0.058286 

-0.074742 

-0. 079605 

C8 

1 

1 

1 

1 

1 

1 

0.055651 

1 

1 

1 

0.028899 

-0.004975 

0.046525 

1^0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

D2 

1. 512805 

1. 458642 

1. 110520 

1. 520587 

1. 175784 

1.495048 

1.051067 

D4 

0. 581078 

0.650948 

0.222962 

0. 704215 

0.260565 

0. 608559 

0.248946 


0.195990 

0.176749 

- - - 

0. 178288 

- - - 

0. 111925 

- - - 

^8 

0. 042675 

0.027054 

1 

L 

0. 018078 

- - - 

1 

- - - 


®-Here d*/h* = 0.05 
^Here d*/h* =0.10 

"^For curvefit A in Figure 2. 5e; Here C^q O.OOO 755 and Cj _2 = -0. 001159- 
'^'^For curvefit B in Figure 2. 5e 
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TABLE III.- DATA POINTS USED IN THE WALL PRESSURE CURVEFITS 



*^Pwall 

0 \ 

4 o 

0 ° 

1.120 

18.8° 

0.918 

57 - 5 ° 

0.281 

56.3° 

- 0 . 394 

75.0° 

-0. 824 

93.7° 

-O.927 

112. 5O 

- 0. 824 

131. 1° 

-0.694 

150.0° 

-0.607 

168. 8° 

-0. 565 

180.0° 

-0. 544 



*"Pwall 

0\ 

175^ 

1 — 1 

175° 

0° 

1.032 

1.202 

1.050 

20° 

0.651 

0. 887 

0. 666 

0 

0 

0.0 

0.095 

-0.025 

0 

0 

-0.575 

-0. 666 

-0. 633 

80° 

-0.726 

-0. 950 

-0.747 

100° 

-0. 596 

-0. 760 

-0. 595 

120° 

-0.476 

-0. 601 

-0. 482 

l40° 

-0.458 

-0. 557 

-0. 4i8 

H 

o\ 

0 

0 

-0. 426 

-0. 558 

-0. 4o6 

180° 

-0.438 

1 

-0. 552 

-0. 592 


T++ 


^Here d*/h* = O.O 5 
^Here d*/h* = 0.10 


'Here d* h* = O.O 5 


"^Data taken from Kawaguti (1955) 

^^Data taken from Grove, Shair, Petersen, 
and Acrivos (1964) 


+++, 


Data taken from Thoman and Szewczyk (1969) 



*"Pwall 


200 

0° 

1.000 

13. 950° 

0.8l4 

20. 925° 

0.642 

27. 900° 

0.451 

34. 875° 

0.232 

4l. 850° 

0. 

48. 825° 

-0.256 

55. 800° 

-0.479 

62. 725° 

-0.721 

69.750° 

-0. 894 

76. 725° 

-0. 991 

85. 700° 

-1.032 

90. 675° 

-1.022 

97. 650° 

-0. 977 

104.625° 

-0. 922 

111. 600° 

-0. 874 

125. 550° 

-O'. 791 

159- 500° 

-0. 731 

153. 450° 

-0.698 

167. 4oo° 

-0.684 

180. 000° 

-0.675 
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TABLE IV.- COEFFICIENTS FOR EQUATION (2-5-5) 


^0 

1.072 

a^ X 10^ 

-8.003 

ai^ X 10^ 

10. 592 

ag X 10^^ 

-6.240 

&Q X 10^^ 

1.752 

a]_Q X 10^*^ 

-0.207 

a^L2 X 10^^ 

0.005 


“^Reported by Dimopoulos and Hanratty (1968) for Re^ = I 74 
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Figure 2.1 Coordinate System 













- - 










• Thoman and Szewczyk (1969) 
See Table III 


Curvefit B 
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Figure 2.5 


Theoretical Displacement 
Bodies 

(a) Re^ = 40 




Figure 2.5 Continued 

(b) Rej = 174 






Figure 2.5 Continued 

(d) Re^ = 177 




Figure 2.5 Concluded 

(e) Re^ = 200 



Figure 2.6 First-Order Displacement Thickness for a Potential Velocity Distribution. 
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Figure 2.8 Surface Shear Stress Distributions (a) Re^j = 40. 
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gure 2.8 Concluded (e) Re^j = 200. 
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Figure 2.12 Comparison of Experimental and Analytical Shear Stress Distributions 
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■essure Gradient Parameter Distribution Using 
at Re^ = 200. 
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Figure 2.18 Inviscid Pressure Gradient Parameter Distribution Using 
Curvefit B at Re^ = 200. 
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Figure 2.19 Inviscid Surface Pressure Gradient Distribution for Curvefit 
B at Re^ = 200. 














Initial Station 



= 12 Flat Plate Case 



Figure 3.3 Infllence of Initial Profiles on Wall Pressure 
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Figure 3, ID Displacement Body Effects on Flat Plate Dispucement Thickness 
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Figure 3.14 Displacement Body Effects on Fut Plate Displacement Thickness 
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I'i = 12.25 Compression Ram> Sruw 
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Figure 3,16 Displacement Body Effects on Wall Pressure - 
No Curvature Corrections 



0.8 


OUTER EDGE OF BOUNDARY LAYER AS 
MEASURED FROM F^OTOGRAPHS 


Y*(lNCHES) 



25 



X 


Figure 3.18 DISPLACE^E^^• Body Effects on Wall Heat Transfer 
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FIGURE 3. 19 EFFECTS OF INVI SC ID CURVATURE CORRECTIONS 
ON WALL PRESSURE 
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H = 12,25 Compression Ramp Stuev 



1st order souurioN 











Figure 3.27 Effect of Free Stream Unit Reynolds Number on Wall Pressures 



X 


Figure 3,28 Effect of Free-Stream Reynolds Number on Heat Transfer to Wall 
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